
Prepared by 

NATIONAL BUREAU OF STANDARDS 
Boulder, Goto, ROJfi 
for Letch Research Center 

JtATiONAl AERONAUTICS AND SPACE AO«!!ffSTtATlON » WASHINGTON/ ftut' • JANUARY W% 


1. Report No. 


2 


. Government Accession No. 


3. 


Recipient's Catalog No. 


4. 


NASA CR-2156 

Title and Subtitle 


5. Report Date 


CAVITATION IN LIQUID CRYOGENS 
H - HYDROFOIL 


December 19 72 

6. Performing Organization Code 


7. Author(s) 


8. Performing Organization Report No. 


J. Hord 

9. Performing Organization Name and Address 

National Bureau of Standards 
Boulder, Colorado 80302 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 

15. Supplementary Notes 

Project Manager, Werner R. Britsch, Fluid System Components Division, NASA Lewis Research 
Center, Cleveland, Ohio 

16. Abstract 

This document constitutes the second of four volumes to be issued on the results of continuing 
cavitation studies. Boundary layer principles, along with two-phase flow concepts, are used 
to improve existing correlative theory for developed cavity data. Details concerning cavity 
instrumentation, data analysis, correlative techniques, and experimental and theoretical 
aspects of a cavitating hydrofoil are given. Both desinent and thermodynamic data, using 
liquid hydrogen and liquid nitrogen, are reported here. The thermodynamic data indicated 
that stable thermodynamic equilibrium exists throughout the vaporous cryogen cavities. The 
improved correlative formulae were used to evaluate these data. A new correlating param- 
eter, based on consideration of mass limiting two -phase flow flux across the cavity interface, 
is proposed. This correlating parameter appears attractive for future correlative and pre- 
dictive applications. Agreement between theory and experiment is discussed, and directions 
for future analysis are suggested. The front half of the cavities, developed on the hydrofoil, 
may be considered as parabolically shaped. 


None 

10. Work Unit No. 

11. Contract or Grant No. 

C-39004B 

13. Type of Report and Period Covered 

Contractor Report 

14. Sponsoring Agency Code 


17. Key Words {Suggested by Author(s)) 

Cavitation Pumps 

C ryogenic s V entur i 

Hydrofoil 
Nucleation 

I 

18. Distribution Statement 

Unclassified - unlimited 

19. Security Classif. (of this report) 
Unclassified 

20. Security Classif. (of this page) 

Unclassified 

21. No. of Pages 
162 

22. Price* 

$3.00 


* For sale by the National Technical Information Service, Springfield, Virginia 22151 













CONTENTS 


Page 

1. Summary 1 

2. Introduction * 3 

3. Analytical Considerations 6 

3. 1 Convection Model (boundary layer analysis) . . 8 

3. 2 Convection Model (dimensional analysis) ... 26 

3. 3 Estimation of V 0 31 

2 cp 

4. Experimental Apparatus 37 

4„ 1 Hydrofoil, Sting -Mount, and Tunnel 37 

4. 1. 1 Design Considerations 38 

4. 1. 2 Details of Fabrication 48 

4. 2 Hydrofoil Contour and Pressure Distribution . 52 

5. Data Analysis 55 

5. 1 Correlation of Desinent Cavitation Data ... 55 

5. 2 Discussion of Desinent Cavitation Data ... 56 

5. 3 Correlation of Developed Cavitation Data ... 66 

5. 4 Discussion of Developed Cavitation Data ... 77 

5. 5 Developed Cavity Shapes 88 

6. Concluding Remarks 89 

7. Nomenclature 91 

8. References 98 

Appendix A: Experimental cavitation data- -nitrogen and 

hydrogen- -for hydrofoil 103 

Appendix B: Thermocouple fabrication procedure 145 

Appendix C: Computer program for correlating desinent 

cavitation data ....... 149 


LIST OF FIGURES 


Page 

Figure 3. 1 -Convective heat transfer model- -thermal boundary- 

layer analysis 10 

Figure 3. 2 Illustration of how choking two-phase flow within 
the liquid- vapor boundary limits mass flow flux 

across the cavity interface 34 

Figure 4. 1 Photographs of hydrofoil and sting -mount 

assembly 39 

Figure 4. 2 Sketch of instrumented, hydrofoil and sting 

assembly 40 

Figure 4. 3 Photographs showing typical appearance of 
vaporous hydrogen cavities on the tapered 

hydrofoil 41 

Figure 4. 4 Details of plastic tunnel 42 

Figure 4. 5 Details of hollow hydrofoil 43 

Figure 4. 6 Details of hydrofoil mandrel 44 

Figure 4. 7 Exploded view of sting assembly 45 

Figure 4. 8 Contour of the cylindrical leading edge of the 

0. 5-caliber hydrofoil 53 

Figure 4. 9 Pressure distribution on the 0. 5-caliber hydrofoil, 

for non-cavitating flow 54 

Figure 5. 1 Effects of tunnel inlet velocity and liquid 

temperature on required inlet pressure for desinent 
cavitation in liquid hydrogen: ( 7 ) = P @ 37 R; 

® = P @ 38. 5R; © = P @ 40R; 

^ v v 

© = P @ 41. 5R 59 

w v 


IV 



LIST OF FIGURES (Continued) 

Page 

Figure 5. 2 Desinent cavitation parameter for liquid hydrogen 

as a function of tunnel inlet velocity and liquid 

temperature 60 

Figure 5. 3 Effects of tunnel inlet velocity and liquid 

temperature on required inlet pressure for desinent 

cavitation in liquid nitrogen: (T) = @ 140R; 

(2) = P @ 1 50R; © = P @ 160R 61 

v v 

Figure 5. 4 Desinent cavitation parameter for liquid nitrogen 
as a function of tunnel inlet velocity and liquid 

temperature 62 

Figure 5. 5 Pressure and temperature depressions within 

cavity in liquid hydrogen 67 

Figure 5. 6 Pressure and temperature depressions within 

cavity in liquid hydrogen 68 

Figure 5. 7 Pressure and temperature depressions within 

cavity in liquid hydrogen 69 

Figure 5. 8 Pressure and temperature depressions within 

cavity in liquid hydrogen 70 

Figure 5. 9 Pressure and temperature depressions within 

cavity in liquid hydrogen 71 

Figure 5. 10 Pressure and temperature depressions within 

cavity in liquid nitrogen 72 

Figure 5. 11 Pressure and temperature depressions within 

cavity in liquid nitrogen 7 3 

Figure 5. 12 Pressure and temperature depressions within 

cavity in liquid nitrogen 74 

Figure 5. 13 Pressure and temperature depressions within 

cavity in liquid nitrogen 7 5 


v 



LIST OF TABLES 


Page 

20 


Table 3. 1 
Table 3. 2 

Table 5. 1 
Table 5. 2 
Table 5. 3 
Table A- 1 a. 
Table A- lb. 
Table A- 2 a. 
Table A-2b. 


Estimated (analytic) exponent values . 

Physical parameters, primary dimensions 
and dimensionless groups resulting from 
dimensional analysis of the cavitation problem . 
Temperature -compensated desinent data 
(Hydrogen: 0. 312-inch ("0. 792-cm) hydrofoil) 

Temperature -compensated desinent data 
(Nitrogen: 0. 312-inch (0. 792-cm) hydrofoil) 
Correlative results for liquid hydrogen and 
liquid nitrogen developed cavity data . 
Experimental cavitation data for 0. 3 12 -inch 
hydrofoil using liquid nitrogen (English Units) 
Experimental cavitation data for 0. 312-inch 
hydrofoil using liquid nitrogen (SI Units) . 
Experimental cavitation data for 0. 3 12 -inch 
hydrofoil using liquid hydrogen (English Units) . 
Experimental cavitation data for 0. 312-inch 
hydrofoil using liquid hydrogen (SI Units) . 


. 28 
. 57 
. 58 
. 78 

. 103 
. 113 
. 123 
. 134 


vx 


CAVITATION IN LIQUID CRYOGENS 
II - HYDROFOIL 
J. Hord 

1. SUMMARY 

This document constitutes the second of four volumes to be 
issued on the results of continuing cavitation studies. The first volume 
dealt with venturi experiments, and this second volume is concerned 
with a 0. 5 -caliber hydrofoil. 

Convective heat transfer and two-phase mass-flow flux limiting 
concepts are used to extend, and improve, existing formulae for cor- 
relating developed cavitation data. These correlative expressions are 
also suitable for predicting the cavitating performance of a particular 
piece of equipment from one fluid to another. Elementary boundary 
layer principles and two-phase flow models are used to derive these 
expressions. Anew correlating parameter, MTWO, is derived; MTWO 
is a liquid phase velocity ratio that is deceptively similar to a Mach 
number and is derived from two-phase flow considerations. 

Details of the hydrofoil-tunnel configuration, cavity instrumenta- 
tion, data analysis, and correlative techniques are discussed. 

Experimental data resulting from this study are presented in 
their entirety in tabular form. Selected data are also presented in 
graphical form. Both desinent and developed cavity data were acquired, 
using liquid hydrogen and liquid nitrogen test fluids. A mathematical 
technique, for correlating the desinent data, was developed and is 
explained in detail. The desinent data exhibit a nearly single- valued 
IL , for each fluid, at the maximum velocities, irrespective of fluid 
temperature. Similar results were obtained with the venturi. The 
desinent hydrogen data were more temper ature -dependent than the 


nitrogen data. Comparison of our desinent data with that of others- -for 

higher boiling-point liquids -- shows that cryogenic liquids require less 

subcooling to avoid cavitation, i. e. , less Net Positive Suction Pressure 

(NPSP) is required for the cryogens. 

Thermodynamic data, consisting of pressure and temperature 

measurements within fully developed hydrogen and nitrogen cavities, 

indicated that stable thermodynamic equilibrium prevails throughout 

the vaporous cavities. These data were correlated using the extended 

theory derived herein. In certain instances, data correlation is 

improved by using liquid kinematic viscosity as a correlating parameter. 

Maximum benefit, in data correlation, is obtained by using the new 

MTWO parameter. When MTWO is used, it appears that only one other 

correlating parameter, cavity length, may be required. Use of MTWO, 

to correlate or predict the performance of cavitating equipment, is 

recommended in all future work. 

Correlative expressions, experimentally evaluated, are in fair 

agreement with those derived from the boundary layer analysis. More 

analysis is needed to explain the subtleties, of the combined mass and 

heat transfer processes, that govern the cavitation process. 

K . is known to vary widely with equipment geometry. This 
c, mm 

study demonstrates that K . can also vary appreciably with flow 

c, mm 

conditions, although the equipment geometry is fixed. This informa- 
tion, coupled with that of others, indicates that current predictive 
techniques must be used with prudence, i. e. , the techniques that rely 

upon constant K . admit an additional source of error for some 
c, mm 

equipment geometries. 

Photographic studies, performed during this experiment, indicate 
that the cavities formed on the hydrofoil have an elliptical shape. The 
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front halves of these cavities are adequately represented by a simple 
algebraic expression of parabolic form. 

2. INTRODUCTION 

Vaporous cavitation is the formation of the vapor phase within a 
flowing liquid, due to a reduction in pressure. Since the formation and 
collapse of vapor cavities alters flow patterns, cavitation may reduce 
the efficiency of pumping machinery [1]^ and reduce the precision of 
flow measuring devices. Collapse of these vapor cavities can also 
cause serious erosion damage [2] to fluid-handling equipment. While 
the noncavitating performance of hydraulic equipment may be predicted 
from established similarity laws, cavitating performance is much more 
difficult to predict from fluid- to -fluid. Recent advances in this area 
have been made by NASA-LeRC personnel [3-6] and others [7-9], but 
additional work is required to improve the current technique for pre- 
dicting cavitating performance of equipment from fluid-to -fluid. The 
effects of fluid properties on cavitation performance are well recognized 
[10-19] and require more understanding to develop improved similarity 
relations [19] for equipment behavior. Much more knowledge is needed 
to extend this predictive capability from one piece of equipment to 
another, i. e. , a more general predictive technique, applicable to equip- 
ment design, must include the effects of equipment geometry and size, 
in addition to fluid properties. 

NASA has undertaken a program [1] to determine the cavitation 
characteristics of various hydrodynamic bodies and the thermodynamic 
behavior of different fluids, in an effort to obtain improved design cri- 
teria to aid in the prediction of cavitating pump performance. The ex- 
perimental and analytical study described herein was conducted in sup- 
port of this program. 


1 

Numbers in brackets indicate references at the end of this paper. 
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Liquid hydrogen and liquid nitrogen were chosen as test fluids for 
this study for the following reasons: (1) the ultimate goal of this pro- 

gram is to acquire sufficient knowledge to permit intelligent design of 
pumps for near-boiling liquids, and (2) predictive analyses indicated 
[1] that the physical properties of hydrogen and nitrogen make them 
particularly desirable test fluids. The objectives of this study were 
1) to experimentally determine the flow and thermodynamic conditions 
required to induce desinent (or incipient) and developed cavitation on 
various hydrodynamic bodies, 2) to improve existing correlative expres- 
sions for the prediction of cavitating performance of hydraulic equip- 
ment, and 3) to establish, if possible, a technique for predicting the 
fluid -handling capability of different cavitating equipment using differ- 
ent fluids. The latter two items are extensions of the state-of-the-art 
and the last objective is highly optimistic, i. e, , accounting for the 
effects of equipment geometry and size in the predictive expressions. 

This report covers the work performed on a two dimensional body 
(hydrofoil). Similar results, using a transparent plastic venturi, were 
presented in Volume I of this report series [20]. 

In the desinent (or incipient) cavity studies, the test section inlet 
velocity was varied from 103, 7 to 200. 5 ft/s (31,6 to 61. 1 m/s) with 
hydrogen, and from 24. 2 to 74. 1 ft/ s (7. 4 to 22. 6 m/ s) with nitrogen. 
Inlet fluid temperatures were varied from 36, 58 to 41. 60 R (20. 32 to 
2 3. 1 1 K) with hydrogen, and from 137.45 to 160.06 R (76. 36 to 88. 92 K) 
with nitrogen. Both incipient and desinent cavitation data were acquired 
with hydrogen and nitrogen test fluids; however, the incipient nitrogen 
data were not used in the final data cor relations -- see appendix C. A 
slight hysteresis effect was detected with nitrogen, but there was no 
noticeable hysteresis in the hydrogen data; i. e. , for hydrogen the flow 
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conditions associated with barely visible cavities did not vary, even 
though the data point was approached from noncavitating (incipient) 
flow, and from developed cavitating (desinent) flow. In this report, 
desinence refers to barely visible cavities, and for convenience the 
terms incipience and desinence are used interchangeably. 

Pressure and temperature profiles, within fully developed cavities, 
were measured and are referred to herein as developed cavitation data. 
Cavity lengths varied from 0. 35 to 1, 75 inches (0 o 89 to 4. 45 cm) with 
hydrogen and from 0. 45 to 2. 00 inches (1. 14 to 5. 08 cm) with nitrogen. 
Test section inlet velocities were varied from 113. 2 to 219. 1 ft/s 
(34. 5 to 66. 8 m/s) with hydrogen and from 26. 7 to 80. 2 ft/s (8. 1 to 
24. 4 m/s) with nitrogen. Inlet liquid temperatures were varied from 
36. 47 to 41. 53 R (20. 26 to 23. 07 K) with hydrogen and from 137. 41 to 
159. 79 R (76. 34 to 88. 7 K) with nitrogen. The bulkstream vapor pres- 
sure exceeds the measured cavity p ressure and the saturation pressure 
corresponding to the measured cavity temperature; therefore, the 
measured pressure depressions within the cavity, and also the pressure 
depressions corresponding to the measured temperature depressions, 
are called "pressure depressions. M Alternatively, the pressure depres- 
sion may be expressed in terms of its equivalent equilibrium "tempera- 
ture depression. M Contrary to the venturi tests [20, 21], no thermo- 
dynamic metastability was detected within the vaporous hydrogen or 
nitrogen cavities developed on the hydrofoil; i. e. , measured tempera- 
tures and pressures within the cavitated regions, appeared to be in 
the r mo dynami c e quilib r ium. 

A similarity equation, based upon the B-factor concept of Stahl 
and Stepanoff [11], has been developed [19] for correlating cavitation 
data for a particular test item from fluid -to -fluid; this correlation is 
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also useful in extending the velocity and temperature range of data for 

any given fluid. This correlative expression has been subjected to some 
2 

criticism, and obviously has some theoretical defects. In this report, 
thermal boundary layer considerations are used to place this correla- 
tive expression on a firmer analytical base; also, the similarity approach 
is extended to include the effects of convective heat transfer. Two addi- 
tional correlating parameters, liquid viscosity and liquid-vapor surface 
tension, emerge from this analysis. A convection-oriented dimensional 
analysis is also performed and a two-phase liquid-vapor sonic velocity 
correlating parameter is introduced. Two-phase mass -flow flux limiting 
concepts are then used to combine this two-phase sonic velocity parame- 
ter, with the liquid inlet velocity parameter. A new and valuable cor- 
relating parameter, MTWO, results. The correlative expressions 
resulting from these analyses are evaluated using experimental data 
from this study. 

3. ANALYTICAL CONSIDERATIONS 
In this section, we shall endeavor to upgrade the theoretical 
framework for correlating developed cavitation data and the cavitating 
performance of equipment. The basic approach of Gelder, et al. [19], 
with some major improvements, will be used; this approach is an 
extension of the B-factor concept introduced by Stahl and Stepanoff [11]. 

A history of the development of the B-factor theory and an improved 
technique for calculating B-factor were recently published [23]; this 
article [23] indicates that the n quasi-static n vaporization model, based 
upon isentropic flashing of the liquid, is reasonably compatible with the 
physical processes of cavitation. Of some concern is the assumption 
that this isentropic vaporization occurs under stable thermodynamic 


See discussions at end of papers by Moore [5] and Ruggeri [22]. 
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equilibrium conditions. There is experimental evidence [20] of slight 
thermodynamic me ta stabilities in the midregions of vaporous hydrogen 
cavities, although most data [19-21] indicate that thermodynamic equili- 
brium exists near the leading edge of vaporous cavities in .a wide variety 
of liquids. Consequently, it appears that the B-factor concept is com- 
patible with the physical processes of cavitation, so long as we restrict 
our attention to the leading edges of the cavity where vaporization most 
surely occurs. The B-factor concept may be viewed as an abstract 
quantity, analogous to entropy in the study of thermodynamics, used to 
obtain a simple solution to our complex and physically real problem. 
Further improvement of the B-factor approach is desirable because 1) 
the concept of isentropic vaporization at the leading edges of the cavity 
is reasonably compatible with the conditions of our problem, 2) this 
approach is well established [19] and provides good results when applied 
to cavitating equipment [4-6], and 3) current efforts to account for the 
dynamic effects of cavitation via bubble growth [9], or entrainment 
theory [24], have not yet resulted in a correlative technique that is less 
complex or less dependent upon experimental data. Also, it appears 
that the B-factor approach consistently provides the best correlative 
results [9]. Thus, it is the intention of this paper to provide a firmer 
theoretical foundation for an improved B-factor correlative approach. 

It will be shown that the simple thermal conduction model used 
by Gelder, et al. , [19] is a special case of a more general 'convection' 
model developed herein. The fundamentals of boundary layer theory 
and convective heat transfer are used to derive the appropriate para- 
metric expressions. Because it is correlating parameters that we 
seek, it is not necessary to concern ourselves with the mathematical 
rigors of an exact and specific boundary layer analysis; however, appli- 
cation of the thermal boundary layer concepts should yield the pertinent 
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heat transfer correlating parameters, thereby placing the B-factor 
approach on firmer ground. Under certain engineering assumptions, 
we shall establish a convective heat transfer model believed applicable 
to the developed cavitation process. In those areas where assumptions 
are required, and little or no information exists, it will be apparent 
that there is a willingness of the author to be guided by experimental 
data. From this analysis, liquid viscosity and surface tension emerge 
as additional correlating parameters. A dimensional analysis confirms 
the boundary layer analysis and introduces the possibility of still an 
additional correlating parameter- -the liquid-vapor two-phase acoustic 
velocity. 

3. 1 Convection Model (boundary layer analysis) 

That convective heat transfer must be involved in the cavitation 

process is apparent. It is therefore equally apparent that the conductive 

i 3 

heat transfer model used by Gelder, et al. , [ 1 9 J is inadequate. Despite 
its shortcomings, the method of Gelder, et al. , has been surprisingly 
successful in predicting the cavitating performance of liquid pumps. 
From this, we can only conclude that their correlative expression, 
though simplified, contains the essential correlating parameters. The 
conduction model used by Gelder, et al. , [19] was based on an analysis 
performed by Eisenberg and Pond [2 5]; this analysis simply relates the 
heat extracted from a liquid film, of constant thickness, to the bulk- 
stream velocity, and to the heat transferred through this film by con- 
duction. Gelder, et al. , [19] also assumed that the cavity thickness did 
not vary with length; it has subsequently been shown [26] that correla- 
tion of experimental data with theory can be materially improved, by 
assuming an appropriate shape for the cavity. In the analysis that 


See discussions at end of paper by Ruggeri [22]. 
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follows, a forced convection heat transfer model is proposed; this 
model features a liquid film with dimensions proportional to the thermal 
boundary layer, and a parabolic-shaped cavity. The resulting expres- 
sions permit the computation of B-factor which is then related through 
the entropy calculations [2 3] to the physical cavitation process, i. e. , to 
the cavity pres sure depression, or temperature depression, and ulti- 
mately to the predicted [3-5] equipment operating conditions. 

The two-dimensional model is shown in figure 3. la; figure 3. lb 
will be used later in considering the axisymmetric case. A number of 
engineering assumptions are required to proceed with the analysis, and 
each is sufficiently important to be discussed in detail. Referring to 
figure 3. la, we will assume that flow is forced over the flat plate by an 
imposed pressure gradient in the bulkstream liquid, i. e, , by external 
means unrelated to the temperature field in the liquid. Body forces 
are hereby neglected, and it will be further assumed that the vapor-filled 
cavity springs from the leading edge of the flat plate (or submerged 
body) --the latter is in good agreement with experimental results. We 
will then assume that the liquid flows over the cavity as if it were a 
solid body. Sufficient heat transfer must occur to vaporize enough 
liquid molecules to sustain the cavity under steady flow conditions; there- 
fore, we will assume the existence of a thermal boundary layer in the 
liquid flowing over the liquid-vapor interface (cavity wall). Because of 
the uncertainty of the flow patterns, circulation, etc. , of the vapor 
within the cavity, it is difficult to assess the virtue of this assumption. 
Certainly, the vaporous cavity must lubricate the liquid flow, and the 
relative velocity of liquid/vapor is unknown; consequently, it is possible 
that the liquid flow at the cavity interface may be laminar even though 
the bulkstream liquid is highly turbulent. It is also possible that the 
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3. la Basic flat plate model. 
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3. lb Axisymmetric model. 


Figure 3. 1 Convective heat transfer model--thermal boundary 
layer analysis. 
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liquid flow at the liquid-vapor interface is turbulent; therefore, we 
shall examine both laminar and turbulent thermal boundary layers in 
the liquid flowing over the liquid-vapor interface. Under these assump- 
tions, the roughness of the interface (cavity wall) can also influence the 
turbulence level and affect the heat transfer coefficient [27 --p. 75]. The 
fluid bulkstream Reynolds number will be turbulent, i. e. , > 10^ in 
most cavitating systems, but the possibility of laminar flow is not to 
be discounted. 

In connection with the foregoing assumption, we shall permit no 
blowing or suction at the interface; in effect, this means that no mass 
transfer from the liquid to vapor phase will be allowed at the interface. 
Obviously, this assumption is violated because the liquid must have a 
velocity component normal to the interface, and there must be some 
mass exchange between the liquid and vapor phases to sustain the cavity; 
however, the liquid/vapor density ratio is large for all liquids, and the 
cavity may easily be sustained by the vaporization of relatively few 
molecules of liquid in the thermal boundary layer. This justification 
is less valid as we approach the fluid critical point and the liquid/ vapor 
density ratio decreases. Again, this f no- suction 1 assumption is difficult 
to evaluate, because we have insufficient knowledge concerning the net 
flow of vapor through the cavity and the rate of vaporization required 
to sustain this vapor flow and hence the developed vaporous cavity. 

With due consideration to the foregoing arguments, and to the other 
assumptions required in this particular analysis, the introduction of 
mass transfer concepts does not appear warranted. This decision 
becomes even more plausible when we recall that we only seek cor- 
relating parameters- -the functional relationship will be established 
empirically. 
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We shall also assume that this limited vaporization, at the inter- 
face, has negligible effect on the shape and thickness of the thermal 
boundary layer and the vaporous cavity. For simplicity, we shall also 
neglect the curvature of the liquid-vapor interface in deriving an expres- 
sion for the thermal boundary layer thickness; this interface curvature 
is retained in describing the cavity thickness. Neglecting interface 
curvature and considering the cavity as a solid body that springs from 
the leading edge of the hydrodynamic body, reduces our problem to that 
of flow over a wedge. The cavity will always originate at or near the 
leading edge of the hydrodynamic body, and the foregoing assumptions 
essentially eliminate the geometrical configuration of the hydrodynamic 
body from our boundary layer analysis. The decision to neglect curva- 
ture of the interface, in the evaluation of boundary layer thickness, does 
not necessarily result in flow over a wedge--we could, for example, 
merely require that the curvature of the interface and the curvature of 
the boundary layer be sufficiently close that transformation to linear 
coordinates does not induce appreciable error. While the shape of the 
hydrodynamic body does not appear to be important to the boundary layer 
analysis, it certainly affects the shape of the attached vapor cavity; 
consequently, the cavity vapor volume and subsequent empirical results 
will be influenced by the body shape. We will therefore restrict our 
correlative results to applications where geometric and dynamic simi- 
larity can be reasonably assured, i. e. , where identical or similar 
hydraulic equipment is used. 

We will assume that the vapor cavity is relatively thin, so that 
we may approximate our problem--for any arbitrary shape--by con- 
sidering forced flow over a flat plate. The thermal boundary layer in 
the liquid will also be taken as sufficiently thin that interaction with 
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other boundary layers can be ignored, i. e. , the body is immersed in 
a semi -infinite flow field (no bounding walls). The hydrodynamic and 
thermal boundary layers will be assumed to originate with the cavity, 
at the leading edge of the body, whether the flow be laminar or turbulent. 
This condition causes us to specify that the liquid-vapor interface be 
at constant temperature [27--p. 204, 234], Experimental evidence 
[19-21] indicates the interface temperature varies slightly, but as 
previously noted we prefer to restrict our analysis to the leading edges 
of the cavity where little temperature variation occurs. 

Consistent with this argument is the assumption that liquid prop- 
erties remain constant as the liquid flows along the cavity wall. Fortu- 
nately, relatively small temperature-depressions, due to vaporization, 
occur along the cavity interface and the fluid properties cannot vary 
significantly. In addition, we have the option of evaluating the fluid 
properties at the reduced (cavity) temperature or at a mixed mean 
temperature if we so desire. Past experience indicates that the physi- 
cal properties can be evaluated at the bulkstream temperature without 
sacrifice in the correlative results. To maintain constant fluid prop- 
erties, we must also neglect liquid compressibility and limit our flow 
to relatively low velocities, i. e. , the effects of kinetic energy and vis- 
cous energy dissipation are neglected. We will require the freestream 
velocity to be invariant along the cavity wall. The hydrodynamic bound- 
ary layer is thicker than the thermal boundary layer in laminar flow, 
but the two are of about the same thickness in turbulent flow T27--p. 234]. 
The latter situation results [27--p. 234] because the transfer mechanism 
through a turbulent boundary layer is primarily eddy diffusion, and the 
eddy diffusivities for heat and momentum are nearly equal where the 
liquid Prandtl number is in the range 0. 5 to 10- -our range of interest. 
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It is assumed that boundary layer separation (stall) does not occur, 
and that the isentropic flashing model [23] is applicable to the vaporiza- 
tion occurring at the leading edge of the cavity. Compatibility of this 
model with the actual cavitation process has already been discussed. 

We would like to emphasize that liquid-vapor thermodynamic equilibrium 
normally exists [19-21] near the cavity leading edge, with metastability 
occurring in some cases [20] in the central regions of the cavity. With 
due respect to this experimental data, we note that radial probing (tem- 
perature and pressure measurement normal to the cavity wall) has not 
been, and perhaps cannot be, accomplished; thus, we have no guarantee 
that stable thermodynamic equilibrium does or does not exist along the 
entire liquid-vapor interface. To the best of our knowledge, it is safe 
to assume that liquid-vapor equilibrium does exist along the leading 
edges of the cavity. Experimental data indicate, as would be expected 
[20], that slower liquid velocities promote liquid-vapor equilibrium 
throughout the developed cavities. 

Under the stipulated assumptions, we cover a wide range of 
practical cavitating flow applications, e. g. , flow over hydrofoils, tur- 
bines, inducers, impellers, cylinders, spheres, ogives, flow inside 
nozzles, etc. For simplicity, we will assume that the 'convective' heat 
transfer to bodies of arbitrary shape may be expressed in the form 

_ _ ^ ___ ml ^ nl , „ , v 

Nu = C Re Pr , (3-1) 

x o x 

where C , ml and nl have different values for laminar and turbulent 
o 

flows. This simple expression is universally accepted and needs no 
further justification. More specific Nusselt expressions, that include 
additional dimensionless numbers, could be used but would merely 
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clutter the analysis presented here; therefore, eq (3-1) will be used to 
complete this analysis. Subsequently, these generalized results may 
be modified to accommodate specific Nusselt expressions --examples 
are given. We will also assume that the characteristic thermal bound- 
ary layer thickness may be expressed in the form 

6 = C x Re“ m2 Pr" n2 , (3-2) 

-tlx ’ 

where again C^, m2, and n2 have different values for laminar and tur- 
bulent flows. 

Using eq (’3-1) and the proposed 'cavity model', it is easily shown 
that eq (3-2) has the appropriate algebraic form for the enthalpy and 
conduction thicknesses of the thermal boundary layer. Later, in this 
section, we shall see that a simplified boundary layer analysis may be 
employed without regard for the numerous foregoing assumptions; how- 
ever, such an analysis does not provide the framework for modifying 
or improving the theory, nor does it shed light on the various physical 
characteristics of the cavitation process. 

Equation (3-2) is the proper result for laminar flow over a flat 
plate, and is consistent with the reasoning advanced by Kays [27--p. 

2 34, 244] for turbulent flow over a body of arbitrary shape. For light 
liquids and small values of n2, the Prandtl term in eq (3-2) approaches 
unity and 6 is approximately equal to the hydrodynamic boundary layer 
thickness for turbulent flow over a flat plate--one of our basic assump- 
tions. Recall that in laminar flow, molecular transport processes 
dominate, while in turbulent flow, eddy transport processes are domi- 
nant. Therefore, it is only in the sublayers of the turbulent boundary 
layer that Prandtl number has any influence. Since the sublayers , 

comprise only a small portion of the turbulent boundary layer thickness. 
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the thermal and velocity boundary layers should have about the same 
thickness regardless of Prandtl number. To maintain generality, irre- 
spective of flow mode, we will retain the Prandtl term in eq (3-2). 

Under the assumptions set forth herein, the generalized solutions 
[27 --p. 222, 244] for heat transfer in laminar and turbulent flow, over 
bodies of arbitrary shape, reduce to the form given in eq (3-1). 

Finally, let us assume that the volume (thickness) of liquid 
involved in the sustenance of the vaporous cavity is proportional to- -not 
necessarily equivalent to--the characteristic thermal boundary layer 
thickness. It seems natural that more liquid can take part in the vapori- 
zation process as particles of liquid progress along the cavity wall, i. e. , 
more time is available for heat to be transferred from liquid located 
further away from the liquid-vapor interface. Thus, we propose mating 
this physical concept of cavitation to the thermal boundary layer develop- 
ment in external flow. We will then require that the liquid-vapor inter- 
face (cavity wall) can be mathematically described by an expression of 

the form 6 = CLx^ . Data from Rouse and McNown \Z8], and data 

v 2 

obtained in this study, indicate that this form is quite adequate and the 
value of p normally lies between 0. 4 and 0. 9 for a wide variety of hydro- 
dynamic body shapes. For this analysis a mean value of p 0. 6 5 will 
be used. Later, we shall use a more general expression for the cavity 
shape to account for body size effects. 

Admittedly, there is no shortage of assumptions to this point, but 
it must be equally obvious to the reader that a rigorous solution to our 
problem is hopelessly complicated. Accordingly, we seek correlating 
parameters, and eventually an empirical solution, for the prediction of 
the cavitating performance of hydrodynamic equipment. Before pro- 
ceeding with the development of the correlating expressions, let us 
tabulate the pertinent assumptions that limit the scope of this analysis: 
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1) Liquid flow is forced over the hydrodynamic body by external 
means unrelated to the fluid temperature field (no body forces). 

2) The liquid flow is steady, incompressible, with relatively low 
and constant free-stream velocity over a hydrodynamic body in a semi- 
infinite flow field. The flow may be laminar or turbulent. 

3) Liquid flow over the vaporous cavity is treated as flow over a 
solid body. The thermal and velocity boundary layers originate with 
the cavity wall, at or near the leading edge of the hydrodynamic body. 

4) The cavity wall is at constant temperature, all fluid properties 
are invariant and are evaluated at the free-stream temperature. Liquid 
Prandtl numbers are <, 10. 

5) Depletion of liquid molecules, due to vaporization, from the liq- 
uid film associated with the sustenance of a vaporous cavity, is assumed 
to have a negligible effect on the shape and thickness of this liquid film 
and the cavity, i. e. , mass transfer across the liquid-vapor interface 

is neglected. Curvature of the liquid-vapor interface is also neglected 
in computing the thermal boundary layer thickness; the interface curva- 
ture is taken as parabolic in computing the cavity thickness. 

6) Boundary layer separation (stall) does not occur. 

7) Equations (3-1) and (3-2) are applicable to flow over developed 
cavities on bodies of arbitrary shape. 

8) The thickness of the liquid film, involved in the sustenance of the 
vaporous cavity, is proportional to the characteristic thermal boundary 
layer thickness. 

9) The "quasi-static" isentropic flashing theory is applicable to the 
cavitation process near the leading edges of the cavity. 

10) The generalized results, obtained from flat plate or wedge-flow 
analyses, are applicable to bodies of arbitrary shape. Date correlation 
is applicable to hydrodynamic bodies with geometric and dynamic simili- 
tude, i. e. , similarity of shapes, fluid motion and forces. 
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The B-factor is defined as the ratio of vapor to liquid volume 
associated with the sustenance of a fixed vaporous cavity in a liquid. 
Then, referring to figure 3. la, we may write 

B =V /V = C X 6 dxA X 6 dx (3-3) 

v t J o V / J o t 

for a plate of unit width (normal to the 6 -x plane). Substituting eq (3-2) 
and the algebraic expression for 6 into eq (3-3), we obtain 


B = [ C 2 fo xPdx ]/[c i Pr' n2 y o X xRe x ' m2 dx] . (3-4) 


Evaluating Re and performing the integration in eq (3-4) produces 


_ _ p — 1 m2 n2 

B = C 0 Re Pr 
3 x 


(3-5) 


For the particular case where m2 = n2, the liquid viscosity is elimi- 
nated from Re and Pr and eq (3-5) reduces to 
x 

B = C.x P " 1 (V x/a) m2 . (3-6) 

3 o 

If we now restrict eq (3-6) to laminar flow, m2 = 0. 5, and if we follow 
Gelder, et al. , [19] and take the cavity thickness as invariant, p = 0; 
performing the indicated algebra on eq (3-6), we get 

B = C,(V /ax)°* 5 . (3-7) 

3 o 


This particular solution is the same result obtained by Gelder, et al. , 
using the conductive model. Had Gelder, et al. , assumed that the 
liquid-vapor interface varied linearly with length (p = 1), they would 
have obtained 


B = C 3 (V q x/a )°* 5 , 


(3-8) 
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a result which has the appropriate arrangement of parameters as 

verified by subsequent experiments [3, 19-21], Note that the product 

of Re and Pr is the Peclet number, Pe = V x/cy, and it is this number 
x * o 

upon which the Gelder, et al. , [19] correlative technique is based. 

Referring back to eq (3-5), we will combine the dimensionless 
numbers to derive the appropriate correlating parameters for the 
general case 


„ _ -n2 m2 p-l+m2 m2-n2 

B = C 3 a V 0 x (p/h) 


(3-9) 


The similarity constant for the cavitating body, , may now be elimi- 
nated and empirical evaluation made possible [19] by relating the pre- 
dicted B to a reference value evaluated from experimental data 



(3-10) 


where El = n2, E2 = m2, E3 = p-l+m2, E4 = m2 -n2 , and v = |i / p , 
the kinematic viscosity. An additional correlating parameter, v, is 
introduced. 

To obtain estimates of the exponent values in eq (3-10), let us 
consider both laminar and turbulent flows and evaluate eq (3-2) in 
three different ways. We will use only the simplest flat plate theory 
and the conventional relations [27] for the conduction, enthalpy, and 
thermal boundary layer thicknesses. Using a mean value for p (^ 0. 65) 
and performing the necessary algebra, we obtain the estimated range 
of values for each exponent as indicated in table 3. 1. 
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Table 3. 1: Estimated (analytic) exponent values. 


Flow 


i2 = El m2 = E2 E3 


E4 Comments 


Laminar 

ml = 0. 50 
nl = 0. 33 
p = 0. 65 


0. 33 
0. 67 


0. 33 


0.50 0.15 0.17 

0.50 0. 32 -0.17 

0.50 0.15 0.17 


Turbulent 

ml = 0, 80 
nl = 0. 60 
p = 0. 65 


0. 60 
0. 40 


0 . 00 


0. 80 


0 . 20 

0 . 20 


0. 45 0.20 

-0.15 -0.20 

-0.15 0.20 


Laminar 


. Turbulent 


0 . 00 


0. 67 


0 . 20 


0. 80 


-0,15 -0,20 Estimated 

to to range of 

0. 45 0. 20 Exponent 

Values 


* - Evaluated from conduction thickness of thermal boundary 

layer. 

* - Evaluated from enthalpy thickness of thermal boundary 

layer. 

* - Evaluated from thermal boundary layer thickness. 


The most physically relevant method of evaluating 6 may be 
indicated by comparison of experimental data with analytical results 
given in table 3. 1. Then, if the proposed model is reasonably valid, 
the heat transfer mechanism that limits heat transfer to the cavity inter- 
face may be revealed. A convection-limited process should be best 


represented by exponents evaluated from the enthalpy thickness, while 
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a conduction- limited process should favor the use of conduction thick- 
ness* The location, within the boundary layer, of the major source of 
resistance to heat transfer will vary with the Reynolds and Prandtl 
numbers. In turbulent flow, where Pr > 1, most of the thermal resist- 
ance resides in the laminar sublayers of the boundary layer; thus, use 
of the conduction thickneoS is suggested. In laminar flow, the thermal 
resistance is distributed over the entire thickness of the thermal bound- 
ary layer, and conduction and convection are of equal importance. 
Because most cavitating flows are turbulent, one is tempted to prefer 
the fourth row of exponents in table 3. 1; however, in reality, the liquid 
flow is vapor-lubricated and the shapes of the hydrodynamic and thermal 
boundary layers are uncertain. Allowing for vapor lubrication- 
irrespective of the liquid flow mode --the liquid velocity must be nonzero 
at the liquid-vapor interface, and the thermal boundary layer must be 
somewhat characteristic of laminar flow* This consideration would 
encourage us to use the enthalpy thickness and expect the second and 
fifth rows of exponents in table 3. 1 to bracket experimental data. 

Obviously, the theoretical values given in table 3. 1 can vary with 
the method of evaluating 6 , cavity shape, and the appropriate Nusselt 
relationship for each specific cavitating system. 

It is also apparent from the foregoing analysis that eq (3-10) 
contains only the most common correlating parameters. Other Nusselt 
relationships, or cavity shape formulations, could easily introduce addi- 
tional correlating parameters as indicated by the following discussion. 

If we were so inclined, we could assume that the heat exchange 
at the liquid-vapor interface retains the salient features of both nucleate 
boiling and forced convection; nucleate boiling is selected because the 
liquid- to -vapor AT is generally within the nucleate regime. Using the 
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superposition theorem. [ 29 ], we would add the individual contributions 
of the boiling and convection processes to obtain the combined heat 
transfer coefficient. This coefficient will decrease, with increasing 
liquid surface tension [29], in a manner dependent upon the contribution 
of the boiling coefficient to the combined heat transfer coefficient. Since 
we are empirically evaluating each correlating parameter in eq (3-10), 
we may introduce surface tension (a) as a simple ratio to test its impor- 
tance as a correlating parameter. If the boiling coefficient strongly 

influences the heat exchange, at the cavity wall, this addition will be 

- E 5 

valuable because the boiling coefficient is ag . Following the develop- 

- E 5 

ment that led us to eq (3-10), we may conclude that B <*o „ Use of 

the boiling coefficient would also introduce other parameters that we 
have not considered, namely vapor density, latent heat of vaporization, 
and appropriate fluid pressures and temperatures. All of these parame- 
ters, excluding a, are in essence used to compute B (hereafter called 
BFLASH) in the isentropic flashing model [23]. To evaluate [20] the 
exponents on eq (3-10), we force B as computed from eq (3-10) to 
match BFLASH as closely as possible. Thus, a is the only parameter 
in the boiling coefficient expression that is unused in the correlative 
expression. Use of the boiling coefficient is not considered entirely 
appropriate because the boiling process differs from the cavitation 
process. The latter is the formation of a vapor phase inside an origi- 
nally jrniiorm_ temperature field, while boiling occurs in an originally 
nonuniform temperature field (surface being hotter than bulk); however, 
a radially nonuniform temperature field surrounds the developed cavity, 
so that the concept of a boiling heat transfer component is tenable. 

Whether the mechanism of vapor generation is mainly attributable to 
a pressure reduction, or to a temperature difference induced by the 
pressure reduction, is of little importance to us. 
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It is also quite likely that o could enter the correlative expression 
through 6 . Dynamic analysis of the stability of the liquid-vapor inter- 
face may reveal a dependency, of the cavity shape, upon the Weber and/or 
Bond numbers. Because the inertia forces are usually much larger than 
the gravitational forces, the Weber number should prevail. The dimen- 
sional analysis performed herein confirms that a enters through the 

2 . 

Weber number. We = p V x/cr is important in the study of bubble for- 

x o 

mation, breakup of liquid jets, liquid surface ripples, etc. , and has 
long been recognized [30] as a potential correlating parameter. Use of 
the Weber number or the boiling heat transfer coefficient introduces 
only one new correlating parameter, a . Equation (3-10) may be written 
to include a as follows 



The numerical value of E5 may be estimated from the boiling heat trans- 
fer data as ^ 0. 5 to 1. 0. 

Gravitational forces may easily be included in this correlative 
expression by permitting the Nusselt number to be a function of the 
Froude number [24]; however, the ratio of dynamic /gravitational forces, 
of the cavitation model, is normally large and this refinement does not 
appear warranted. Note that use of the Bond number in describing the 
shape of the cavity would introduce both gravity and a. 

Next, we direct our attention to the body size effects that were 
neglected in assuming that 6 = C x^ . Under our assumption of similar 

V L* 

dimensions, velocities, forces, etc. , we can follow the logic of Moore 
and Ruggeri [31] and assert that cavity thickness and cavity length will 
vary directly with scale size, i. e„ , with some characteristic dimension, 
D, of the cavitating body. Then, similar cavitation performance occurs 
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with similar bodies of different size when x/D is held constant, Moore 


and Ruggeri [31] show that this ’model similarity’ reasoning merely 

E6 

introduces a (D/D J term into eq (3-11). In developing eq (3-5), 

rei 

we assumed that 6 = adequately represents the shape of the 


v 


cavity wall, irrespective of body shape or size. We also assumed that 
the cavity wall originates near the leading edge- -actually near the mini- 
mum pressure point--of the hydrodynamic body; then it is obvious that 
the radius of curvature of the cavity wall, near the point of origin, 
varies with the radius of curvature of the body. Geometrical analysis 

of a family of par abolic- shaped curves indicates that we can satisfy 

E G 

this size requirement by replacing C with C / D , where E6 (like p) 

Ld L* 

may have a particular value for each body shape. Since we are neglect- 
ing the effect of cavity wall curvature when evaluating the thermal bound- 
ary layer thickness, the size effect will not enter into the liquid volume 

E6 

calculation, and BaD --see eq (3-4) and eq (3-5). Then, accounting for 
size, eq (3-11) may be revised 

B = B 


ref 




(3-12) 


The size term is identical to that obtained by Moore and Ruggeri [31]. 

The foregoing examples clearly indicate that the constituent 
parameters of the chosen Nusselt expression, and cavity shape formula- 
tion, must appear in the correlative expression, eq (3-10), 

So far, our analysis has been restricted to two-dimensional flow. 
Referring to figure 3. lb, we consider the axisymmetric body. Rewriting 
eq (3-3) for the axisymmetric body produces 



Rearranging terms, we obtain. 



6 + 2R 
v 



6+26 + 2R 
l v 



(3-14) 


Comparing eq (3-14) with eq (3-3), we note that some second-order 
terms have appeared; therefore, for first order approximations, eq 
(3-14) reduces to eq (3-3) 0 Also, recalling that our problem deals 
with the cavity wall near its leading edges, we may argue that in most 
cavitating equipment 2 R » 6 and 2 R » 6 ; thus, eq (3-14) reduces 

V 

to eq (3-3)„ In consideration of the various assumptions used in this 
analysis, the refinement offered by eq (3-14) does not appear warranted; 
so we shall assume that eq (3-3) and subsequent expressions, through 
eq (3-12), are applicable to cavitating flow over plane-surface and 
axisymmetric bodies* 

Before turning to dimensional analysis, it should be pointed out 
that the equations developed in this section may be derived in a simpli- 
fied manner* It is only necessary to assume that the heat transfer at 
the cavity interface can be expressed in the general form given by eq 


(3-1), irrespective of body shape* No other assumptions concerning 

the ’cavity model’ are necessary. Next, the conduction thickness of 

the thermal boundary layer is given [27--p. 50] as k/fi ; then 6 <*k / h 

x x 

Evaluating eq (3-1) to obtain h , and combining h with the foregoing 
relation for 6 , results in a functional expression for 6 that is identi- 
cal to eq (3-2) 0 Referring to eq (3-4), we note thatj*6 dx does not con- 
tain any physical parameters^ that are not contained in 1/ J*6^dx. 

Then, to simply obtain correlating parameters, we may state that 
B cc (1/k )J*fi dx 0 So long as we are empirically evaluating exponents 


With the possible exception of D, for scale effects, and D can be 
introduced as per reference [3l] 0 
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on each physical parameter, i. e. , the initial values of ml, nl , m2, 
n2, etc. , are unimportant; this particular expression will produce the 
same results as the generalized approach outlined herein. 

From a practical viewpoint, B<*(l/k) J* ft dx is all that we need 
to develop correlating parameters, but no other information is gained 
from this simplified approach. From a theoretical viewpoint, the gen- 
eral approach outlined herein provides detailed information concerning 
the 'cavity model'. The latter approach, after comparison of theory 
and experiment, may be modified to account for theoretical deficiencies-- 
no such avenue exists for the simplified case. In the simpler case, only 
empirical adjustment of eq (3-1) is possible, and no physical explana- 
tion for the adjustment is visible; however, a variance in one or more 
of the numerous assumptions used in the detailed analysis may well 
synthesize experiment and theory. Although the detailed approach is 
preferable in a theoretical sense, the simplified view that B (l/k)j*ft dx 
will be used to good advantage in our dimensional analysis. 

3. 2 Convection Model (dimensional analysis) 


The primary strength of dimensional analysis is the formulation 
of dimensionless groups of variables, useful in the correlation of experi- 
mental data. A complete analytic solution does not result, but signifi- 
cant dimensionless numbers are identified, e. g. , Reynolds number, 
Prandtl number, etc. The main weakness of dimensional analysis is 
that no information is gained on the inner mechanisms of the physical 
problem--thus, the functional dependency of each dimensionless number 
must be determined empirically. To apply dimensional analysis, it is 
necessary to know which parameters influence the problem, and the 
success of the method depends upon the intuitive selection of these 
parameters. One’s confidence in the results is usually bolstered if 


26 



there exists a skeletal theory and/or some understanding of the physical 
problem to support the parameter selection. Because of the work in the 
previous section and existing correlating techniques [5, 7-9, 19, 22], we 
are in a good position to perform such an analysis. To determine the 
number of independent dimensionless groups required to describe the 
problem, we use the Buckingham tt - theorem and the ML 8 T dimen- 
sional system--M = mass, L = length, 9 = time and T = temperature. 

The Buckingham method produces n groups, tt , tt tt . . . tt , of 

1 2 3 n 

dimensionless variables and predicts a solution to the problem has the 
form tt^ — f (tt^, . nj; n is determined by subtracting the number 

of primary dimensions, e. g. , M, 6, T, from the total number of phy- 
sical variables. Table 3. 2 summarizes the physical variables, pri- 
mary dimensions and resultant rr - groups for the cavitation problem. 

The physical variables listed in this table are either geometric 
parameters or else they may be determined from appropriate pressure 
and temperature measurements. The first seven parameters listed are 
the fundamental convection variables. The eighth parameter introduces 
size effect, and the ninth parameter permits surface tension to influence 
the stability and heat exchange phenomena at the liquid-vapor interface. 
The tenth and last parameter is included because apparent 'shocks' 
have been experimentally observed [20] in cavitated regions of liquid 
hydrogen, and such 'shocks' are known to affect heat transfer [27--p. 

271]. Also, at least two authors [7, 8] have developed correlation 

5 

techniques based upon two-phase compressibility phenomenon. Note 
that the first seven physical variables in table 3. 2 require only four 


5 

See reference [32] also. 
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Table 3. 2 Physical parameters, primary dimensions and 

dimensionless groups resulting from dimensional 
analysis of the cavitation problem. 


Physical 

Primary 

Dimensionless Groups of Variables 

Parameter 

Dimensions 

Resulting from Dimensional Analyses 





X 

L 

1 


V 

o 

L/e 



k 

ML/e 3 T 

f Nu 

1 X 




Nu 

P 

M/L 

\ 

x 

C P,£ 

2 2 
l /e t 

( Re 

1 X 

Re 

X 

M* 

m/l e 

1 Pr 

Pr 

h 

3 

m/ e t 

/ 

x/D 

X 



D 

L 


We 

X 

a 

M/e 2 


M Zcp 

v_ 

Zcp 

L/e 

J 





primary dimensions-- then. 

we anticipate that three tt - groups will be 


required to correlate the data. Buckingham's rule assumes that 


TT = 


a b c d e f 
x V k p C (jl 

o p, 



(3-15) 


where the exponents are to be evaluated to determine the individual 
tt- groups, TTj, tt^, etc. Substituting the primary dimensions of each 
physical variable into eq (3-15), and then equating the sum of the expo- 
nents of each primary dimension to zero, produces a set of equations. 
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For tt to be dimensionless, the sum of the exponents on each primary 
dimension must be zero. This procedure produces four equations, with 
seven unknowns; therefore, three of the exponents may be evaluated 
arbitrarily, so long as they are linearly independent. Simultaneous 
solution of these equations, in a fashion that assures the exponents are 
linearly independent, provides one of the tt - groups. This process is 
repeated twice to evaluate tt and rr s Table 3. 2 shows that this method 

u J 

produces tt = Nu , tt = Re , and rr = Pr; from experience we know 
1 x 2 x 3 

that Nu = f (Re , Pr). 
x x 

If we include size, surface tension, and two-phase sonic velocity 
in the list of physical variables and repeat the dimensional analysis, 
we get 

Nu = f(Re , Pr, x/D, We , ). (3-16) 

x V x x 2 cp 

The result is dimensionless groupings of our input parameters. With- 
out prior theoretical work, we would have no idea of the functional 
relationship of the individual n - groups; however, to conform with the 
Nusselt relationship, we will assume that 

Nu = (Re ) a (Pr) b (x/D) C (We ) d (M_ f , (3-17) 

xx x 2cp 

where the exponents are to be evaluated empirically. By evaluating 
each tt- group in terms of physical parameters, rearranging, and com- 
bining exponents on each physical parameter in eq (3-17), we obtain 

h = ku" b V t x U v V (p/a) d (1/D) C (1/V_ ) 6 , (3-18) 

x o 2 cp 

where t = a + 2d + e, u = a+ c + d-l and v = a-b. Recalling that 
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B cc(l /k) dx and proceeding as in the development of eq (3-10), we 
obtain 



Because p is contained in both the a and v terms, there is little need 

to preserve p as an explicit parameter in eq (3-19). Alternatively, we 

could retain k/C and p terms and eliminate cy. To minimize the 
p, l 

number of correlating parameters and thus simplify the correlative 
expression, we elect to retain a , a proven correlating parameter. The 
resultant expression is very similar to eq (3-12) , 



( 3 - 20 ) 


where the exponents in eq (3-19) have been converted to those used in 
eq (3-12) and are to be evaluated empirically. So, without further 
struggle, the basic correlating parameters are produced, although 
some of the parameters are inverted. 

The term emerges because of our assumption that it is a 

significant physical variable. This parameter may be related to the 
’’characteristic velocity” discussed by Holl and Wislicenus [30] and 
thought to be important in cavitation studies. Jakobsen [8] and Spraker 
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[7] have enjoyed some success in the correlation of inducer and pump 
cavitation data using two-phase acoustical parameters; however. 

Chive rs^ is skeptical of the role of the acoustical parameter in the 
correlation of cavitation data. In this paper, we include the acoustical 
parameter in the B-factor expression, but do not rely upon it as a sole 
or primary correlating parameter. Computation of requires deter- 

mination of fluid quality, a quantity that varies from zero to 100 percent 
across the liquid-vapor interface. Therefore, evaluation of V is not 
straightforward, and the effect of V on and ultimately on B is not 
clear. A method of estimating is proposed in the following section 

of this paper. 

The dimensionless groups in eq (3-17) may ultimately simplify 
correlation of cavitation data, but more insight into the mechanism of 
the problem is obtained by correlating on the individual physical vari- 
ables. Thus, until further knowledge evolves, it is suggested that 
eq (3-20), or a modification thereof, be used for correlative purposes. 
Evaluation of exponents in the correlative expression is fully explained 
elsewhere [20], and use of the results for predicting cavitating per- 
formance of equipment has also been published [3-6]. 

3. 3 Estimation of V 

2cp 

Computation of can only be performed in an abstract manner. 

This situation occurs because varies with fluid quality and the 

assumed thermodynamic process. For our problem, we cannot define 
the fluid quality or the thermodynamic process. Of the various thermo- 
dynamic models available [33], it is proposed that the homogeneous 


6 


See discussion at end of papers by Chivers [9]. 
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metastable model is the most appropriate. This model is selected 
because it works best [33, 34] for experimental systems with fluid 
qualities < 0. 10--our range of interest, as will be demonstrated. Also, 
it is based upon an isentropic flashing process, which is consistent with 
the BFLASH theory [2 3], In this model, the liquid is assumed incom- 
pressible and an equivalent two-phase pressure -volume relationship, 
for an isentropic process, is used to derive the following relation [33], 


Z 


G 

c 


ZC + C 

EiZ XjA 

ZC + c 

v, v v, t 


P P 

1 v ( Z+l ) 

z 


(3-21) 


where Z = X/(l-X). Using the basic definitions of X and B, we readily 
obtain the algebraic expression 


X = 


B Z_ 

B + (p ,/p ) = Z + 1 ’ 

U V 


(3-22) 


and Z = Bp /p 
v l 


( 3-Z3) 


Then eq (3-Z1) may be written as a function of fluid quality, X, 

or BFLASH, as we see fit* For our problem we have no independent 

way of estimating X but a proven method of estimating BFLASH [23]. 

G (= V ) is the choking (limiting) mass flow flux for two -phase 
c 4cp Zcp 

flow, based upon the fluid quality as computed from eq (3-ZZ). Thus, 
it is apparent that the two-phase acoustic velocity, Y , a f unc ti° n 
of BFLASH. Now, since BFLASH is a characteristic vapor/liquid vol- 
ume ratio, for a prescribed cavity, it is clear that G is a parameter 

c 

that is characteristic of this cavity and cannot vary locally along the 

cavity interface. Therefore, G ( = p V ) may be used to compute a 

c 

characteristic liquid velocity component, V , normal to the cavity 

L 

interface. This liquid velocity component may be regarded as the 
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characteristic limiting rate at which liquid may approach the liquid- 

vapor interface. This velocity component is not to be confused with the 

actual limiting rate at which liquid crosses the interface, but should be 

considered as representative of that actual rate. Then, in this simplified 

and abstract fashion, we can account for mass transfer considerations 

that were neglected in the foregoing analyses. Figure 3. 2 illustrates 

this mass flow flux limiting concept, i. e. , mass transport to the liquid- 

vapor interface is limited by choking two-phase flow across the interface. 

Combining the characteristic mass -limiting liquid velocity, V , with the 

f, 

liquid inlet velocity, V , forms the velocity ratio, V /V . This velocity 

o of, 

ratio will be designated MTWO, to emphasize that this ratio is similar 

to a Mach number and is proportional to the two-phase Mach number, 

M^ , i. e. , MTWO = V /V and M^ = V /VT , then MTWO = M^ V /V . 
2cp o f, Zcp o 2cp 2cp 2cp f. 

Because G = p_ V\ = p o V . and VT /V_ = p „/p 0 , we may alternatively 
c 2cp 2cp l V 2cp l l 2cp 7 7 

state that MTWO = M_ p o / p_ . Also, V 0 = G / p o and 

2 cp f 2 cp f, c f, 

MTWO = V /V =V p /G . (3-24) 

of, o f, c 



Reasons for evaluating (and correlating on) the MTWO parameter, 
rather than the (or parameter, will be discussed in detail in 

a later section of this report. 

Several thermodynamic models are available [33] for computing 
G and a fundamental relationship for acoustic velocity in a homogeneous 
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Figure 3. Z Illustration of how choking two-phase flow within the 

liquid- vapor boundary limits mass flow flux across the 
cavity interface. 
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compressible fluid was derived by Jakobsen [8], Using the Jakobsen 
formula and eq (3-24), we obtain 


V 


MT WO = 


1 + Z(p /p ) 2 (a /a ) 2 

t V l V 

1 + Z 


0. 5 


(3-26) 


Equation (3-25) is restricted to relatively low fluid qualities 
( < 0. 10), but eq (3-26) bears no such restrictions,, These two equations, 
evaluated with experimental data, produce almost identical trends (not 
equivalent values); therefore, they produce almost identical correlative 
results when used with the correlative expression, eq (3-20). Because 
eq (3-26) bears no fluid quality limitations and is no more complex than 
eq (3-25), it is recommended for future use. Introducing MTWO into the 
correlative expression of eq (3-20), and eliminating the V and V 

O Ct cp 

terms, we obtain 


B = B 


ref 



MTWO 


MTWO . 

ref 



(*)>)» 


E6 


(3-27) 


Subsequently, we shall show that this expression is far superior to 
eq (3-12) in correlating the experimental cavitation data acquired in 
this study. 

The two-phase thermodynamic models advocated by Smith [33], 
which are the vapor- choking, the homogeneous thermal equilibrium, and 
the homogeneous metastable models, have been tested with experimental 
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data from this study. Of these models, only the homogeneous metastable 
model produced beneficial correlative results. This model, evaluated 
with formulations for both compressible and incompressible fluids, cor- 
related our experimental data with equal success. To test the validity 
of the two-phase-limited mass transport concept, the experimental data 

were correlated 1) using V /V , 2) using V /V and 3) using V /a . 

o l oZcp o v 

Beneficial correlative results were obtained only when V /V was used, 

o l 

substantiating the choice of MTWO as a correlating parameter. Note 
that eq (3-25) and eq (3-26) are evaluated at the reduced pressure and 
temperature conditions at the cavity interface - -in this study, these 
reduced fluid conditions are measured. BFLASH is also evaluated at 
the reduced fluid temperature. 

It will be convenient, in a later discussion, to express eq (3-26) 

in a simplified form. This is accomplished by first noting that 

p / p » B for all of our experimental data, i. e. , B < 5 and p / p > 40. 
I v ^ t v ^ 

Then, from eq (3-22), X ^ B p /p and from eq (3-23), X % Z. Also, 

v Kj 

X ^ Z < 0. 1 3 for our experimental data--the choice of the homogeneous 

metastable two-phase model is made apparent. Referring to eq (3-26), 

/ 2 3 

we note that (p^a^/ p^a^) is a very large quantity ( > 10 ); therefore, 

for values of Z of interest (^0.01 to 0. 13), we may simplify eq (3-26) 

as follows 
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V « (1/BFLASH) 

As 


(3-30) 


0. 5 

These two relations will be useful in explaining data correlations, 
in another section of this report. 

4. EXPERIMENTAL APPARATUS 
The experimental apparatus used in this study was explained in 
detail in the first volume [20] of this report series. The experimental 
facility, instrumentation, error statements, visual and photographic 
aids, and test procedures are fully described in that document [20]. 

One additional error statement-- concerning uncertainty in pressure 
measurement for nitrogen test fluid- -is needed and given in section 5. 2 
of this report. Only the test section (tunnel) and hydrodynamic body 
(hydrofoil) details need to be discussed here. This tunnel was located 
between the supply and receiver dewars of a blowdown flow system, 
see reference [20]. 

4. 1 Hydrofoil, Sting-Mount, and Tunnel 
The 0. 5-caliber hydrofoil, used in this experiment, was chosen 
so that developed cavitation test data could be obtained for external flow 
over a two-dimensional body. These data may ultimately be correlated 
with similar data for external flow over axisymmetric bodies (ogives) 
and for internal flow through a venturi. Also, the hydrofoil offers the 
ultimate opportunity to study the effect of cavity shape on the correlative 
formulae, i. e. , the effects of flow as symmetries, fluid rotation, etc. , 
should be minimized and the cavity thickness, as a function of length, 
velocity, etc. , should be easily determined. 

The two-dimensional body (hydrofoil) used in this experiment 
was a tapered plate with a cylindrical leading edge (0. 5 caliber). Ini- 
tially, the hydrofoil was designed as a flat plate with a cylindrical 
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leading edge; however, difficulties were encountered with tunnel cavi- 
tation, and it was necessary to slightly taper the plate. Using the maxi- 
mum plate thickness, the length-to-thicknes s ratio of the hydrofoil is 
8:1 . A transparent plastic tunnel, and a rigid metallic sting-mount, were 
designed to experimentally complement the 0* 5 -caliber hydrofoil. Photo- 
graphs of the hydrofoil and sting-mount assembly are shown in figure 4. 1. 
A sketch of the instrumented hydrofoil and sting assembly is given on 
figure 4. 2, and figure 4. 3 shows photographs of the hydrofoil as installed 
and viewed during a test. Details concerning the tunnel, hydrofoil, and 
sting -mount are given below and on figures 4. 4 through 4. 7. 

4. 1. 1 Design Considerations 

The hydrofoil was initially designed as a flat plate with sufficient 
tunnel blockage to cause the hydrofoil to cavitate readily; however, pre- 
liminary tests indicated that cavities, developed on the hydrofoil, were 
sufficiently thick to cause tunnel cavitation. The tunnel (test section) 
cavitation is attributed to the pressure reduction that accompanies the 
acceleration of liquid between the cavity and tunnel walls. Selection of 
a thinner hydrofoil alleviates this Bernoulli effect in the vicinity of the 
cavity, but simultaneously decreases the hydrofoil cavitation number, 

K. --resulting m -. lower magnitudes of subcooling, P - P , in the tunnel 

IV O V 

inlet liquid, and therefore enhancing cavitation in the inlet and down- 
stream rep : -ns of the tunnel. Thus, to perform experiments within the 
pressure an .low limitations of the existing facility, it was necessary 
to taper the hydrofoil, f ront-to-rear. This 1. 817 degree taper provided 
sufficient relief to avoid tunnel cavitation, without detachment of the 
cavity on the hydrofoil, i. e. , the cavity developed on the hydrofoil 
adhered to the taper--all tests were monitored via remote closed- 
circuit TV and motion picture cameras. 
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Fi-Sure 4. 1 Photographs of hydrofoil and sting-mount assembly. 
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Figure 4. 3 Photographs showing typical appearance of vaporous 
hydrogen cavities on the tapered hydrofoil. 
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Figure 4 # 5 Details of hollow hydrofoil, 




OGIVAL END-CAP 



Figure 4. 7 Exploded view of sting assembly. 


The slight taper required to eliminate tunnel cavitation could 
have been provided by tapering the sidewalls of the tunnel; however, 
this approach was rejected because of the difficulties in machining a 
smooth, non-wavy, divergent passage (diffuser) in a square cross 
section acrylic tunnel. This plastic tunnel was designed with a square 
flow passage to enhance two-dimensional flow and to promote smooth 
transition to circular cross sections--inlet and outlet ducts were cylin- 
drical, This square tunnel permits direct observation of cavity shape 
and thickness with minimum optical distortion- - a primary objective 
of this study. Consistent with this objective was the necessity to design 
a snug fit between the edges of the hydrofoil and the tunnel, so that 
cavitation would not occur prematurely at the hydrofoil-tunnel interface. 
This was easily accomplished by taking advantage of the differential 
thermal contraction between the acrylic tunnel and the brass hydrofoil. 

A snug fit at room temperature assures a tight fit at cryogenic tempera- 
tures; however, the fit must not be sufficiently tight to fracture the plas- 
tic tunnel when cooled to cryogenic temperatures. The square corners 
of the tunnel flow passage were rounded, to reduce stress concentra- 
tions, and to discourage cavitation along these corners. This rounding 
facilitates cleaning and polishing of the tunnel. The square -to -circular 
transitions, at the inlet and outlet of the tunnel, were designed to elimi- 
nate cavitation at these locations. Also, the inlet transition was designed 
to promote uniform two-dimensional flow at the hydrofoil. Uniformity 
of flow in the tunnel inlet was experimentally verified by static pressure 
measurements at different locations in the tunnel inlet, see figure 4. 4. 

The sting-mount assembly was designed to 1) eliminate cavita- 
tion in the vicinity of the sting -mount and 2) provide a rigid support for 
the hydrofoil (and subsequently for scaled ogives). The latter criterion 
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was satisfied by providing a close -tolerance tongue-and-groove joint 
between tunnel and hydrofoil, see figures 4. 4 and 4. 5. This tongue- 
and-groove also provided assurance that the hydrofoil was properly 
positioned in the tunnel. The hydrofoil-sting-tunnel assembly was 
designed so that the hydrofoil-tunnel interface completely shields the 
tongue-and-groove joint, thus avoiding premature cavitation along this 
joint. The aft (downstream) end of the sting was fitted with a hollow 
4. 33-caliber ogive, see figure 4. 7. This ogive was attached by threads 
to permit ready access to the instrument wires and tubes, inside the 
cylindrical sting, and to facilitate installation of this instrumentation. 
This ogive -shaped cap was hollowed out to allow room for thermally 
anchoring the thermocouple wires. 

The thermocouple junctions and pres sure -sensing ports, used to 
measure temperatures and pressures within the cavities developed on 
the hydrofoil, were located as shown in figure 4. 5. All of the sensors 
were located on one side of the hydrofoil- -with pressure and tempera- 
ture sensors at each station closely spaced- -to ensure a fair compari- 
son of the pressure profiles obtained from the pressure and tempera- 
ture measurements within the cavity. Although the cavities developed 
on the hydrofoil (top and bottom sides) appeared symmetrical, this 
approach assures that both the pressure and temperature sensors are 
subjected to the same cavity environment. Also, the sensors were 
located as near as possible to the central plane of the hydrofoil, to 
avoid potential 'edge effects,' within the cavity, at the tunnel-hydrofoil 
interface. 

The hydrofoil test assembly was installed in the same space 
allocated for the plastic venturi [20] in the experimental apparatus; 
fabrication details follow. 
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4, 1. 2 Details of Fabrication 

The tunnel was constructed from annealed, cast acrylic rod. 

A 31/32-inch (2.46 cm) diameter hole was drilled through the axis of 
this cylindrical rod and reamed to 63/64-inch (2. 50 cm). This hole 
served as a pilot for a square-shaped broaching tool; thus, the square- 
hole passage was broached using a three-pass broaching operation. 

The first broaching tool rough-cut the square hole. The second broach- 
ing tool finish- cut the square hole, and the third broaching tool burnished 
the passage. The broaching tools were specially designed, by an indus- 
trial tool manufacturing firm, to broach plastic with minimum waviness 
and maximum smoothness. Broaching an acceptable square hole required 
a special technique. It was necessary to complete each broaching opera- 
tion with a single thrust, i. e, , the broaching tool had to be shoved all of 
the way through the plastic without stopping. The rate at which the broach 
was shoved through the plastic was not critical. Satisfactory broaching 
operations were performed using a vertical drill press, hydraulic press, 
and an ordinary machine lathe. The only requirements were that the 
machine provide smooth, continuous feed, and sufficient travel to com- 
plete the broaching operatipn in one stroke. Several tool lubricants were 
tried (e. g. kerosene and light machine oil), and we learned that a com- 
mercial additive for automotive engine oils worked best. 

Following the broaching operation, the square tunnel was slotted, 
on two opposing sides, to provide the groove for the tongue -and-groove 
hydrofoil- tunnel joint. This slotting was accomplished with an anchor- 
shaped tool, and an axial scraping motion, on a milling machine. The 
square-to-circular transitions, on the ends of the tunnel, were then 
made by making a taper -cut with a machine lathe. This machining opera- 
tion leaves tapered gullies at each corner of the square passage, see 
figure 4.4. The edges of these tapered gullies were then very carefully 
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rounded, by hand, using appropriate Swiss pattern files and polishing 
compounds. The entire internal passage was then carefully polished 
to a high lustre, using plastic polishing compound. To ensure uniform 
polishing on each sidewall, a square-shaped wooden bar was machined. 
This bar, when covered with a thin cotton cloth, was designed to extend 
through, and fit comfortably inside of, the tunnel. Polishing compound 
was applied to the cloth, and polishing was accomplished by grasping 
the wooden bar on each end and sliding it back and forth in the tunnel. 

With the critical machining and polishing completed, the tunnel 
exterior was machined and polished. Internal tunnel dimensions were 
easily measured by conventional means. Tunnel dimensions are given 
on figure 4. 4. The locations of pressure taps, upstream and down- 
stream of the hydrofoil leading edge, are also shown. Note that the 
exterior surfaces of the tunnel walls are machined flat and parallel 
to the interior surfaces, so that cavities may be observed with mini- 
mum optical distortion. Scribe marks on the tunnel were used to esti- 
mate developed cavity lengths, see figure 4. 3. 

The sting assembly consisted of a hollow coaxial cylinder, 
supported by four radial struts, see figure 4. 7. Two of these struts 
were hollow ellipsoids, and two were very thin solid blade- struts with 
rounded edges. The hollow struts served as conduits, for pressure 
transmitting tubes and thermocouple wires, as did the hollow cylindrical 
sting. These sting mount features are shown on figures 4. 1 , 4. 2, and 
4. 7. The sting was soldered into place, in the mounting flange, as 
shown on figures 4. 1 and 4. 2. Coaxial alignment of the sting was 
maintained, through the soldering operation, by alignment pins and 
the slot-inserted blade-struts, see figures 4. 1 and 4. 2. Alignment 
of the hydrofoil, when soldered to the sting socket, was guaranteed 
by using a specially devised assembly fixture. 
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Five pressure transmitting tubes were routed from the hydro- 
foil, through the hollow sting, and up through one of the hollow ellip- 
tical struts. Then, the tubes pass through a slot in the flange and extend 
through a short length of 0. 25-inch (0. 64 cm) diameter tube that is 
located outside of the flange. The smaller tubes are collectively soldered 
inside this larger tube, to form a seal, and then the larger tube is 
sealed to the flange with a commercial compression fitting. Similar 
fittings are provided for the pressure sensing stations on the plastic 
tunnel. The small pressure transmitting tubes terminate, outside of 
the flange, in the vacuum insulation space. These tubes are then solder- 
connected to larger tubes, which penetrate the vacuum barrier, and are 
attached to pressure transducers. 

The other hollow elliptical strut is used to guide thermocouple 
wires into the vacuum space. Details of the thermocouple fabrication, 
seals, etc. , are given in appendix B. 

Construction of the hydrofoil is quite intricate; therefore, one 
must carefully study figures 4. 2 and 4. 5-4. 7 to fully appreciate the 
finer details. The hollow hydrofoil is easily manufactured by conven- 
tional means. The 1. 817 degree taper was easily attained with a sine- 
plate attachment for a milling machine, and the cylindrical leading edge 
was obtained by scraping the hydrofoil with a concave -radius wheel- 
cutter; the hydrofoil was held in a vise, and the wheel-cutter was 
mounted in a milling head and drawn across the hydrofoil- -without 
rotation of the cutter. 

The transition from cylindrical sting to hydrofoil is accomplished 
with two diametrically opposed fins and a conically tapered cylindrical 


Stainless steel tubes, 0. 040-inch (0. 102 cm) diameter, with 
0. 005-inch (0. 013 cm) wall thickness. 
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shaft extension, see figure 4. 7. The conically tapered hollow shaft is 
soldered onto the male portion (mandrel) of the two-piece hydrofoil. 

The fins are extensions of the tongue machined into the sides of the 
hydrofoil, and help to gradually decelerate the liquid at the aft end of 
the hydrofoil. 

The mandrel portion of the hydrofoil was designed to accommo- 
date ten of the 0. 040-inch (0. 102 cm) tubes, see figure 4. 6. The ten 
tubes were installed in five tube tracks --two tubes to the track, i. e. , 
one tube on top of another. Each tube was plugged with silver solder 
on the end nearest the leading edge of the hydrofoil. The top tube is 
slightly shorter than the bottom tube--after the hydrofoil is assembled, 
sensor ports are drilled through the hydrofoil, into each tube, penetrat- 
ing only one tube. With the tubes in place, the mandrel and tubes were 
tinned with soft solder. The hollow hydrofoil, as shown on figure 4. 5, 
was puddled full of soft solder. The mandrel and hollow hydrofoil 
were then carefully assembled, while the solder was liquid. 

Upon completion of this assembly, holes were drilled through the 
hydrofoil into each of the ten small tubes as indicated on figure 4. 5. 

Care must be taken to avoid drilling through the top tube into the tube 
which lies below. Each of the holes, so drilled, are isolated from the 
others by the soft solder that fills all voids between the mandrel and 
hollow hydrofoil. Five of these holes become pressure sensing stations, 
while the other five are used for thermocouples. This entire assembly 
is then attached, by soldering, to the sting socket using a special align- 
ment fixture, see figure 4. 2, The small tubes must be threaded through 
the sting, hollow strut, etc. , to mate the hydrofoil and sting. Also, the 
five tubes used as thermocouple conduits are extended through a tube 
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support, at the rear of the sting, and soft soldered to this support. 
Later, during thermocouple fabrication, this support is epoxied to the 
sting. 

The thermocouples were installed as described in appendix B; 
following this installation the thermocouple sensors extend through 
small mounds of epoxy. These epoxy bumps must be removed so that 
1) the bare thermocouple junction is flush with the surface of the hydro- 
foil and 2) the thermocouple junction is surrounded by epoxy that elec- 
trically and thermally isolates the thermocouple while sealing it to the 
hydrofoil. This installation technique assures rapid response of the 
thermocouples, while electrically and thermally isolating the junctions 
from the metal hydrofoil. The epoxy bumps were finished flush with 
the hydrofoil, by using Swiss files and then fine-grit sandpaper; during 
this hand-finishing operation, the entire hydrofoil was tape -masked, 
exposing only the epoxy bumps. Following this finishing operation, 
the entire hydrofoil was carefully polished, measured, and installed 
in the test apparatus. 

4. 2 Hydrofoil Contour and Pressure Distribution 

The actual and theoretical contours of the 0. 5 -caliber hydrofoil 
are shown on figure 4. 8. The actual contour was verified by using an 
optical comparator (shadowgraph), with 31.25 X magnification. The 
theoretical noncavitating pressure profiles for this tapered hydrofoil, 
and a non-tapered hydrofoil, are shown on figure 4. 9; experimental 
data from this study are also plotted for comparison. The theoretical 
pressure profiles were computed using existing computer programs 
[35, 36]. Figure 4. 9 indicates good agreement between experimental 
and calculated data, in the regions where experimental data are avail- 
able. Also, note that there is little difference in the calculated data 
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for the flat and tapered hydrofoils. The slight wiggle, in the cal- 
culated C for the tapered hydrofoil, occurs at the flat-to -diffuser 
P 

transition. Figure 4.. 9 indicates that the minimum pressure point 
occurs at about 78 degrees of arc, measured from the stagnation point 
of the hydrofoil. This minimum pressure location is indicated on the 
contour shown on figure 4. 8. The value of 78 degrees corresponds 
favorably with the pressure coefficient data for a cylindrical body 
with a hemispherical head [28]. 

5. DATA ANALYSIS 

The desinent (incipient) and developed cavitation data, for liquid 
hydrogen and liquid nitrogen, are given in complete detail in appendix 
A. These tabulated data are given in English and metric units. No 
attempt has been made to separate the desinent and developed cavity 
data; however, the desinent (or incipient) cavity data are clearly 
marked by the attachment of asterisks to the run numbers. The 
desinent and developed cavity data are correlated, and discussed 
separately in this section. 

5. 1 Correlation of Desinent Cavitation Data 
With the blow-down facility used in this experimental study, it 
was impossible to maintain a constant fluid temperature, while varying 
the inlet velocities and pressures to obtain desinent cavities. Conse- 
quently, it was necessary to develop a mathematical technique for cor- 
relating the desinent data. This was accomplished by using a least- 
squares surface -fitting computer program. Once an equation is 

obtained, to fit the experimental surface (P , V , T coordinates), 

o o o 

the conventional isotherm data for desinence are readily calculated. 
Complete details concerning this correlating technique, and the com- 
puter program, are given in appendix C, Also, the polynomial 
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expressions, derived to correlate the desinent hydrogen and nitrogen 
data, are given in appendix C. These expressions were used to com- 
pute the desinent data presented in tables 5. 1 and 5. 2. These same 
data are plotted on figures 5. 1 to 5. 4. 

The experimental values of K. are tabulated under the heading, 

iv 

K , in appendix A, along with the experimental values of K for 
developed cavities. The K. parameter, as given in appendix A, in 
tables 5. 1 and 5. 2, and on figures 5. 2 and 5. 4, is not corrected for 
tunnel blockage. Multiplying by the square of the appropriate 
area ratio, A^/A^, corrects for tunnel blockage, i. e. , 

, 2 

Corrected K. = (Experimental K. ) (A /A^) , (5-1) 

xv iv 1 2 

where A^ = blocked cross sectional flow area and A^ = unblocked 

(inlet) cross sectional flow area. This correction factor evolves from 

v 

consideration of steady volumetric flow, C , and Bernoulli’s equation, 

P 

and is derived in detail elsewhere [37, 38]. This correction factor 
assures that the minimum static pressure is the same for blocked 
and unblocked flows, when the freestream static pressures are iden- 
tical. Thus, tunnel constraint is easily accounted for, so that the 
results of this study may be readily compared with other experimental 
data. In this study, the correction factor has a numerical value of 
0.47. 

5. 2 Discussion of Desinent Cavitation Data 

Figure 5. 1 indicates that P varies almost linearly with V , 

o o 

at constant T , for hydrogen test fluid; also, P increases with T , 
o o o 

at constant V , in a conventional manner. The incipient (desinent) 
cavitation parameter, K , for hydrogen also behaves in a conven- 
tional manner, see figure 5. 2, On both of these figures, the 
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Table 5. 1. Temperature-compensated desinent data 
(Hydrogen: 0. 312-inch (0. 792-cm) hydrofoil) 
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Table 5. 2. Temperature-compensated desinent data 
(Nitrogen: 0. 312-inch (0. 792-cm) hydrofoil) 



00 

oo 

00 

oo 

oo oo 

CO 

CO 

CO CO 

CO 

co 

O' 

o 

O' O' O' 

O' 



r- 

r- 

I s - r- r- 

CO 

CO 

CO CO 

CO 

CO 

00 

00 

00 

00 

oo 

00 


• 

• 

• 

• 

• • 

• 

• 

• • 

• 

• 

• 

• 

• 

• 

• 

• 

r- 

I s * 

r- 

r- 

h- r*- 

CO 

CO 

CO CO 

CO 

CO 

00 

05 

oo 

oo 

oo 

00 



r- 

I s - 

r- r- 

oo 

oo 

00 00 

00 

oo 

00 

oo 

00 

00 

00 

00 


_ N 

xO 

7 

7 

8 

8 

9 

vO 

7 

7 

8 

oo 

O' 

6 

7 

7 

00 00 

O' 

rn 

• 

• • 

• • • 

• 

• • • 

• 

• 

• * • 

• • 

• 

5 

r- 

O CO 

rH rH 

16 

19 

22 

I s - 

O 00 sO 
H rH rH 

61 

22 

7 

10 

13 

16 

19 

22 


1 

00 

CO O' 

<r 

o o 

m 

oo 

rH <n 

xO 

O' 

CM CM 

CM 

CM CM 

CM 

£ 

cn 

O CM 

rH 

x0 v0 

iH 

rH 

00 o 

00 

CM 

x* 00 

00 

xO 

4- 

( J 

• 

• • 

• 

• • 

• 

• 

• • 

• 

• 

• • 

• 

• • 

• 


xO 

rH 00 

00 

o m 

r- 

rH 

r- p- 

00 

on 

O' CM 

oo 

r- oo 

CM 


f — 1 

CM CM 

CO 

m vO 

CM 

cn 

cn 

in 

r- 

cn 


m xO 

oo 


V0 



h CO fO fl H 

CM Cn -4* LTk vO 

• • • • • 

CM CM CM CM CM 


rOsOHCMoOsj 
vf xO LT\ LO in \0 

• • • • • • 

cn CNJ CM CM CM CM 


O' I s " 00 CO O) 
C co M to <f in 

• •«••• 

CO CM CM CM CM CM 


O 

ft 


nj 

• H 

co 

a* 


vO H CO (\| O CO 
r- IT\ o CO CO CM 
• ••••• 

co oh in co id 
cm co <r m r- & 


00 (M CO CM r- o 

CO CM CO CM CO CO 
• ••••« 

o m >j- oo in v0 
co in \£> oo o 


oooor^cM^- 
h h oo (\i in in 

cm o co o O' 
m so o oo O' r-t 




O 

o 

o 

0 

0 

o 

o 

o 

0 

0 

o o 

o 

0 

0 

0 

o 

o 

CO 

• 

• 

• 

• • 

m 

• 

• 

• • 

• • 

• 

• • • 

• 

• 


in 

in 

m 

in in 

in 

in 

in 

in in 

in m 

in 

in in m 

n 

n 

-H 

SH 

CM 

cn 

<r 

m xo 

r- 

CM 

m 

xt in 

xO 

CM 

cn m 

xO 

n- 



o 

o 

o 

o o 

c 

O 

o 

o o 

o c 

o 

o o o 

o 

o 


o 

o 

o 

o o 

o 

O 

o 

o o 

o o 

o 

o o o 

o 

o 


• 

• 

# 

• • 

• 

• 

• 

• • 

• • 

• 

• • • 

• 

• 

o 

o 

o 

c o 

o 

o 

o 

o o 

o o 

o 

o o o 

o 

o 


<*■ 

>* 

<r 

<!■ 


in 

n 

in in 

in in 

xO 

xO xO x0 

xO 

xO 


rH 

rH 

rH 

rH rH 

rH 

rH 

rH 

rH rH 

rH rH 

rH 

rH rH rH 

rH 

rH 


58 




v 0 .ft/s 


Effects of tunnel inlet velocity and liquid temperature on 
required inlet pressure for desinent cavitation in liquid 
hydrogen: <T) = P v @ 37R; © = P @ 38. 5R; © = P @ 
40R,<4; = P @41. 5R. v 


Figure 5. 1 
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Figure 5. Z Desinent cavitation parameter for liquid hydrogen as 

a function of tunnel inlet velocity and liquid temperature. 
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Figure 5, 4 Desinent cavitation parameter for liquid nitrogen as a 
function of tunnel inlet velocity and liquid temperature. 
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boundaries of the experimental data are indicated by the solid lines, 

i. e. , the broken-line curves are extrapolations of the experimental 

data. Mathematical extrapolation, by correlative expressions, of 

experimentally-determined surfaces is rather risky; therefore, the 

broken-line extrapolations must be used with great caution. Only 

within the boundaries of the experimental data are the mathematically 

derived data--as presented in tables 5. 1 and 5. 2 and on figures 5. 1 

to 5. 4- -considered totally valid. 

Similar comments may be made about the nitrogen test data 

shown on figures 5. 3 and 5. 4; however, the P , Y , T curves-- 

o o o 

figure 5. 3--are non-linear and figure 5.4 requires additional discus- 
sion. Although the P , V , T curves, on figure 5. 3, show the appro- 

o o o ' ~ 

priate temperature dependency, the K. curves, on figure 5. 4, do not. 

This inconsistency in the curves is believed to be mainly due to 

instrument error, even though this error is very small. In reference 

[20], it was estimated that the maximum experimental uncertainty 

2 

in P , for hydrogen, is ± 1. 0 psi (± 0. 69 N/cm ); for nitrogen, the 
maximum uncertainty of the pressure measurements, including cali- 
bration and read-out errors, is estimated to be within ± 1. 5 psi 

2 2 
(± 1. 04 N/cm ), with an allowance of ± 0. 3 psi (± 0. 21 N/cm ) for 

2 

systematic error and ±1.2 psi (± 0.8 3 N/cm ) for imprecision. For 

very low values of V (and P ), the maximum uncertainty in P , for 

o ° 2 ° 

hydrogen, did not exceed ± 0. 5 psi (± 0. 35 N/cm ). At similar con- 
ditions, for nitrogen, the maximum uncertainty in P did not exceed 

2 

±1.0 psi (± 0. 69 N/cm ). Then, the pressure measurement error, 
for the nitrogen data, is considerably larger than that of the hydrogen 
data- -twice as large at the minimum values of V and P for each 
fluid, and fifty percent larger at the maximum values of and P^ 


63 



for each fluid. The instrument error, for the hydrogen data, was 
smaller because the tests were performed at lower pressures, 
permitting the use of lower range pressure transducers--these trans- 
ducers (and data recording system) have an uncertainty of about 
± 0. 5 percent of full range, irrespective of the range used. 

In reference [21], it was shown that the parameter can 
easily magnify instrument error. This amplification of instrument 
error results because K. is computed from a small pressure differ- 

IV 

ence (P - P ), between two relatively large numbers. Because P 
o v v 

is constant, for a given T , the difference (P - P ) is minimum at 

o o v 

the lowest V , see figure 5. 3. Consequently, the maximum error in 

K. will always occur at the minimum velocity, for any specified fluid 
xv 

temperature. Thus, the unusual behavior of the K curves, on figure 

5. 4, at the lower velocities, is not too surprising. To emphasize this 

point, we extract the following data for run number 294J, as tabulated 

in appendix A- la: T ^ 140 R, V jsa 30 ft/s and P - P ^ 12 psi. Now, 

o o o v 

if we assume that there is no error in p , V , or P , and that the error 

o o v 

in P is ± 1. 0 psi, K. will be in error by ^8. 3 percent; therefore, the 

o ^ IV 

140 R isotherm, on figure 5. 4, could easily be in error by 8. 3 percent 
at the minimum velocity. Obviously, such a displacement could drasti- 
cally alter the shape of the curve. Similar computations, at the maxi- 
mum velocities, where P - P is maximum, indicate that instrument 

o v 

error should be almost negligible. Then, we must conclude that the 

K curves, at the maximum velocities, are very near their true posi- 
iv 

tions; thus, it appears that the K. , for nitrogen, is only weakly tem- 
perature dependent, and approaches the value of 2. 6 at the maxi- 
mum velocities, irrespective of fluid temperature. The unusual shapes, 
and positions, of the K. curves on figure 5. 4 are attributed to this weak 
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temperature dependence, amplification of instrument error, and 

mathematical correlation of the experimental data. Because of the 

uncertainties of the extrapolated P , V , T data, at low velocities, 

o o o 

on figure 5. 3, this data is not plotted on figure 5. 4. The data shown 
on figure 5. 4 are considered of real significance only at velocities 
greater than ~ 45 ft/s (~ 14 m/s). 

Similar difficulties were not encountered with the hydrogen data, 
because of a stronger temperature dependence and lower instrument 
error. The strong temperature dependence, of hydrogen, is apparent 
in figure 5. 2, even though the temperature variation is small in com- 
parison with the nitrogen data shown on figure 5. 4. 

The hydrofoil desinent cavity data are quite similar to the data 
obtained for the quarter-round venturi [21]. For both bodies, the 
nitrogen and hydrogen data tend toward a single -valued K , for each 
fluid, at the maximum velocities. 

The hydrofoil data, presented herein, may also be compared 
with data for a similar hydrofoil (untapered), tested in water [30]. To 
make this comparison, the K data on figures 5. 2 and 5. 4 must be 
multiplied by 0. 47 to correct for tunnel blockage. Direct comparison 
of these data indicates that the asymptotic K. for water is greater 

IV 

than the asymptotic K. for either nitrogen or hydrogen (by about 20 

IV 

to 2 5 percent). Similar results were obtained in the venturi studies 

[21], where the data for nitrogen and hydrogen were appreciably 

lower than the K. data for refrigerant-114 [39], ethylene glycol [40], 

and water [40]. Because the latter three liquids have higher bulk 

densities than the cryogens, it is apparent that the cryogens require 

less subcooling for desinent cavitation to occur, i. e. , for K to be 

smaller for the cryogens, at the same V , requires smaller values of 

P - P --see the definition of K. on figure 5. 4. 
o v IV 
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5. 3 Correlation of Developed Cavitation Data 

Typical profiles of measured pressure depression, for liquid 
hydrogen, are given on figures 5. 5 to 5. 9; similar profiles, for liquid 
nitrogen, are plotted on figures 5. 10 to 5. 13. From these figures, it 
can be observed that, within data accuracy, stable thermodynamic 
equilibrium exists throughout the cavity. This topic is discussed in 
considerable detail in section 5, 4 of this report. 

The correlative expressions, given by equations (3-12) and (3-27), 
were used to correlate the developed cavity data from this experiment. 
These expressions, along with the correlative technique developed by 
Gelder, et al. [19], the isentropic BFLASH theory [23], and two least- 
squares data-fitting computer programs [20], were used to correlate 
these hydrofoil data. For convenience, we will refer to these correla- 
tive expressions as 'similarity 1 equations. The similarity equations 
are used to correlate developed cavitation data, in similar test items, 
and to predict the cavitation performance of a test item from fluid-to- 
fluid, and from one temperature to another, when limited test data 
from a single fluid are available. 

Complete and detailed descriptions of the correlative technique, 
computational steps, and computer programs are given in reference 
[20], The correlative procedure, as previously described [20], can 
be followed directly when using eq (3-12). To use eq (3-27), simply 
substitute MTWO for in the computer program. Briefly, this cor- 
relative procedure ensures that the B values calculated from eq (3-12), 
or eq (3-27), and the BFLASH values [2 3] for each data point, are as 
nearly identical as possible; because both B values, at each data point, 
are evaluated from experimental data, this correlative procedure 
produces the best possible agreement between experiment, the isen- 
tropic flashing theory [23], and the correlative expres sion- -eq (3-12) 
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Figure 5. 5 Pressure and temperature depressions within cavity in 
liquid hydrogen. 
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Figure 5. 6 Pressure and temperature depressions within cavity in 
liquid hydrogen. 
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Figure 5. 7 Pressure and temperature depressions within cavity in 
liquid hydrogen. 
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Figure 5. 8 


Pressure and temperature depressions within cavity in 
liquid hydrogen. 
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Figure 5. 9 Pressure and temperature depressions within cavity in 
liquid hydrogen. 
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Figure 5. 10 Pressure and temperature depressions within cavity in 
liquid nitrogen. 
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Figure 5. 11 Pressure and temperature depressions within cavity in 
liquid nitrogen. 
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Figure 5. 12 Pressure and temperature depressions within cavity in 
liquid nitrogen. 
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Figure 5. 1 3 Pressure and temperature depressions within cavity in 
liquid nitrogen. 
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or eq (3-27). This ’best-fit, 1 of the experimental data, is obtained by 
selecting appropriate exponents for each of the correlative parameters 
in the correlative expression--eq (3-12) or eq (3-27). The exponent 
selecting process is quite complex and is treated in appropriate detail 
in reference [20]. Exponents for eq (3-12) and eq (3-27) were derived 
to evaluate the suitability of MTWO as a correlating parameter. 
Attempts to correlate the experimental data, using eq (3-20), were 
abortive, as will be explained in section 5. 4 of this report. 

In eq (3-12) and eq (3-27), the cavity lengths were evaluated at 
the visually observed lengths. BFLASH was obtained, for each experi- 
mental data point, using the average measured cavity pressure depres- 
sion (P - PJ, T , and the calculation method outlined in reference 
v 1 o __ 

[23], The average minimum cavity pressure, P^ = < P 1 + P 1,T )/2 ' 
was used because the temperature measurements are considered just 
as accurate as the pressure measurements. 

In eq (3-12) and eq (3-27), the fluid physical properties are 

evaluated at P and T , with the exception that MTWO is evaluated at 
o o 

the minimum measured cavity pressure, P^ . The standard deviation 
in B is computed for each set of exponents; the individual exponents 
may be held constant or chosen by the computer. The standard devia- 
tion in B factor is minimized in the computer programs when one or 
more of the exponents is selected by the computer; the absolute mini- 
mum standard deviation is obtained when all of the exponents are 
selected by the computer--as in this report. The standard deviation 
is simply computed in those cases where the exponents are held con- 
stant. The set of exponents that produces minimum standard deviation 
in B is selected as the best correlative solution for any particular 
batch of data; i. e. , the standard deviation is a measure of the validity 
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of the similarity and isentropic-flashing theories, as both are evaluated 
from experimental data. Because MTWO proved to be a valuable cor- 
relating parameter for the hydrofoil data, and it had not been used to 
correlate the venturi data [20], both the hydrofoil and venturi data were 
separately correlated using the approach described herein. The results 
of this effort are given in table 5, 3 and are discussed in the following 
section of this report, 

5. 4 Discussion of Developed Cavitation Data 

Photographs of fully developed vaporous cavities, in liquid hydro- 
gen, are shown on figure 4, 3. Inlet velocity and liquid temperature 
were observed to have very little effect on the appearance of cavitating 
hydrogen; i. e 0 , the cavities were generally well defined and uniformly 
developed. Only at the highest temperature, and lowest velocities, did 
the hydrogen cavities exhibit a slightly porous, non-uniform character. 
Similar, though more pronounced, features were observed in the nitro- 
gen cavities. If a cavity is sufficiently porous, i. e. , if it visually 
resembles vapor streams, erratic developed cavity data may result. 

This occurs, because the pressure and temperature sensing ports are 
not continuously covered with vapor during a test, but may be inter- 
mittently covered with vapor and then liquid. This effect is believed 
to be negligible for all of the hydrofoil data acquired and reported 
herein. 

In figures 5. 5 to 5. 13, the data points representing cavity pres- 
sure measurements have been connected with a smooth curve to facili- 
tate comparison with the data points obtained from the cavity tempera- 
ture measurements. The pressure depressions obtained from the cavity 
temperature measurements are, for the most part, in good agreement 
with those derived from the measured pressures, i. e. , within the 
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allowances of instrument error, the cavity vapor is in stable thermo- 
dynamic equilibrium c Evidence of metastable vapor was presented in 
the venturi study [20]* In figures 5. 5 to 5. 13, the only significant 
differences, in the measured pressure and temperature profiles, occur 
near the trailing edges of the cavities. Similar characteristics were 
observed in the venturi studies. The trailing edges of a cavity are 
normally irregular and are characterized by randomly- spaced clouds 
of condensing vapor. Hence, it is inadvisable to attempt to interpret 
the cavity data in this ill-defined region; however, the following dis- 
cussion will shed some light on this topic. 

The data on figures 5. 5 to 5. 13 indicate that some of the cavities 
were shorter than their visual (as observed on film) length. Because 
the pressure depression should be near zero at the trailing edge of the 
cavity, the actual length of the cavity may be estimated by extrapolating 
the data to zero pressure depression. The actual length of the cavity, 
and the observed length, differ because of the irregular trailing edges 
of the cavity and the difficulty in judging the visual length. Both actual 
and visual cavity lengths were used to correlate the data, and they 
produced essentially the same results. The visual cavity lengths, as 
tabulated in appendix A, were used in the final correlative data fits 
reported in table 5. 3. 

Now, referring to figure 5. 8, we will attempt to explain the 
apparent discrepancy, in pressure and temperature measurements, 
near the aft end of the cavity. Note that both the actual and visual 
cavity lengths indicate that the T^ and T^ sensor stations are down- 
stream of the cavity; however, T and T may well be in a region where 

4 5" 

small clouds of vapor are continuously or intermittently condensing. 

If these temperature sensors were exposed to liquid only, downstream 
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of the cavity, we would anticipate that they would indicate near-zero 

pressure depression, i. e, , the liquid temperature should be very near 

T , so that P - P __ is very small. Then figure 5. 8, and T and T , 
o v n, T 4 o 

indicate exactly what we would expect--a gradual transition from aft-end 
cavity vapor temperature, T^, to bulkstream liquid temperature, 
(represented by zero pressure depression on figure 5. 8). The presence 
of vapor, in any quantity, would tend to cool these downstream tempera- 
ture sensors to a temperature intermediate to T_ and T . Similar 
comments apply to the other pressure depression plots. Referring to 
figures 5 # 5 to 5. 13, we note that the agreement, between the pres- 
sure and temperature measurements, is exceptionally good in those 
cases where the cavity covered all sensors, e. g. , see figures 5. 5, 5. 7, 

5. 10, and 5. 11. Plots similar to figures 5. 5 to 5. 13, of all the data, 
reveal that stable thermodynamic equilibrium exists within the cavity 
vapor--within data accuracy. 

’Similarity* equations, used to correlate cavitation performance 
of a particular flow device from fluid -to -fluid, were fitted with numeri- 
cal exponents derived from the experimental data (hydrofoil and venturi) 
of this study. These equations were derived herein and represent exten- 
sions of the work of Gelder, et al. [19]. The numerical exponents for 
hydrogen-refrigerant 114 test fluids in a venturi were obtained from 
the literature [3]; these exponents are compared, in table 5. 3, with 
those obtained in our venturi studies [20], and with those deduced from 
our hydrofoil experiments. The exponents given in table 5. 3 were 
obtained with a least- square s fitting technique and a digital computer; 
the suitability of the various exponents to the experimental data is 
indicated by the standard deviation in B-factor as explained previously. 

In the hydrogen venturi study, the value of B ranges from two to five; 
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in the hydrofoil experiments, the value of B ranges from one to five, 
for both hydrogen and nitrogen test fluids. 

The expressions used to correlate the experimental data are 
given at the bottom of table 5. 3. As indicated in table 5. 3, there is 
little difference in the exponents obtained for the venturi data, using 
different test fluids, i. e. , the exponents in lines 4 and 5 are in good 
agreement, excluding the diffusivity term. The lack of variation in 
a (< 10 percent) explains why El tends to a negative number (line 4) 
when correlating with liquid hydrogen alone. Gelder, et al. [19], 
obtained a value of 0. 5 for El when correlating refrigerant 114 data 
with 25 percent variation in o\ There was over 400 percent change 
in O' in the hydrogen-refrigerant 114 data correlated by Moore and 
Ruggeri [3], and thus the value for El reported in line 5 of table 5. 3 
is to be preferred when correlating with eq (3-12). The mathematical 
technique, used to derive the exponents, can easily pick an extraneous 
value for any of the exponents if there does not exist significant varia- 
tion in the corresponding physical parameter. 

Use of the viscosity and surface tension terms in eq (3-12) did 
not significantly improve the venturi hydrogen data correlation, and 
consequently exponents for these terms are not included in table 5. 3; 
however, we shall demonstrate that viscosity is an important corre- 
lating parameter when attempting to correlate data from fluid -to -fluid. 
It is anticipated that both of these parameters may improve data cor- 
relation where sufficient variation in the fluid properties occurs. In 
the venturi hydrogen data, the viscosity varied less than 20 percent 
and the surface tension varied less than 40 percent. 

Exponents on the hydrofoil data, using eq (3-12), are given in 
the first three lines of table 5. 3. Again, the viscosity and surface 
tension terms had little influence on the data correlation for single 
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fluids (lines 1 and 2) and were not used; however, the viscosity term 
significantly improved the combined fluid correlation (line 3), and the 
corresponding exponent, E4, was determined. For the combined fluids 
(line 3), inclusion of the viscosity term reduced the standard deviation 
by 1 5 percent--a substantial reduction. In all of our venturi and hydro- 
foil data, use of the v and a terms improved the correlation; however, 
it is felt that use of these additional correlating parameters is not 
justified, unless they substantially improve the correlative fit. There- 
fore, it is apparent that only the combined fluid hydrofoil data 'line 3) 
benefited from the use of v » Correlation of the hydrogen- ref rige rant 
114 data (line 5) would most likely be improved by using one or both 
of these terms. In the nitrogen data (line 2), a varied by only 12 percent, 
V varied by 35 percent, and a varied by 30 percent. In the hydrogen 
data (line 1) a varied by 12 percent, v varied by 16 percent, and a 
varied by 20 percent. Thus, it is not surprising that the v and a terms 
were of little benefit in the single fluid correlations, nor that the expo- 
nent on the o: term seeks a negative value. In the combined fluid data 
(line 3) a varied by almost 100 percent, v varied by 35 percent, and 
a varied by 300 percent. Then the a exponent, El, should be quite 
meaningful, the v exponent, E4, (though beneficial) is suspect, and the 
a exponent, E5, should be beneficial. Because the a term was of negli- 
gible value in the correlative fit, even though it varied by a factor of 
three, we must conclude that a is not an important correlating parame- 
ter for the cryogens tested; however, it may yet prove to be a valuable 
correlating parameter for other fluid combinations - -with smaller or 
larger variations in a. Although the \j term improved the combined 
fluid correlation, the numerical value of E4 is suspect because of the 
relatively small variation in v for these data. Again, \j may be an 
excellent correlating parameter for other fluid combinations and is of 
considerable value for the hydrogen-nitrogen combination. 
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It is interesting to note that the combined fluid correlations (lines 

3 and 5) produced exponents, El, EZ, and E3, of nearly the same value c 
Certainly, it is not expected that the exponents, for different bodies, 
should be the same. Recent tests [24], on zero-caliber ogives, indicate 
that these exponents may vary widely from one hydrodynamic body to 
another. 

Comparison of the best experimental exponents (lines 3 and 5 of 
table 5. 3), with the exponents predicted from heat transfer considera- 
tions (table 3, 1), is not too gratifying. It appears that the predicted 
values of E2 and E3 bracketed the experimental data, but the predicted 
values of El and E4 did not. These tests do not reveal that a particular 
flow mode, and technique for evaluating the thermal boundary layer 
thickness, are to be preferred; however, lines 2 and 4 of table 3. 1 are 
in the best general agreement with the experimental data (lines 3 and 
5 of table 5. 3), Actually, this agreement between experiment and heat 
transfer theory is not too bad, when we consider the numerous assump- 
tions involved in the analysis --particularly the neglect of mass trans- 
fer across the liquid-vapor interface. The following discussion indi- 
cates that mass transfer must be considered. 

The venturi and hydrofoil data were also correlated using eq 
(3-27); this was accomplished by substituting MTWO for in the com- 
puter program. These results are shown on lines 6 through 9 in table 
5. 3. It is apparent from the foregoing discussion, and the results 
shown in table 5. 3, that the v and a terms did not materially improve 
the correlation. That eq (3-27) is quite superior to eq (3-12), as a 
correlative expression, is readily shown by direct comparison of the 
results given on lines 1 and 6, lines 2 and 7, lines 3 and 8, and lines 

4 and 9* A marked reduction (~ 1/3) in standard deviation of B-factor 

is achieved, in each case, by substituting MTWO for V • A line-by-line 
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comparison of these same data will also show that the importance of 
the a, and x terms is generally lessened when eq (3-27) is used, i. e # , 
the numerical values of El and E3 are reduced. The numerical value 
of El is sometimes so small (lines 8 and 9) that the a term could, in 
these cases, be neglected. The cv term is undoubtedly diluted some- 
what, because of the many thermophysical and thermodynamic fluid 
properties embodied in the MTWO parameter. It is believed that the 
slight correlative improvement offered by the v and ct terms, when 
using eq (3-27), can be attributed to the very strong influence of the 
MTWO parameter. While the v term was of considerable benefit (line 
3), when using eq (3-12), it has little effect when the standard deviation 
is reduced to a much lower value (line 8), by use of the MTWO parame- 
ter. The predominant influence of the MTWO term may also be respon- 
sible for the slight reductions in E3, when using eq (3-27). The impor- 
tance of MTWO, as a correlating parameter, emphasizes that mass 
transfer plays an important role in the cavitation process. 

An attempt was made to correlate the experimental data using 
both V and (or ) as se P ara ^e correlating parameters, i. e, , an 
expression of the form given in eq (3-20) was used. This attempt 
failed because V °c (1/BFLASH)°* 5 , see eq (3-30). Now B FLASH will 
appear on the right-hand side of eq (3-20) because V is a function of 
BFLASH; also, the exponents on all of the terms are computed by mini- 
mizing the differences in calculated B (left side of eq (3-20)) and BFLASH 
at each data point. Then the computer merely selects an exponent, for 
the V term, that causes B (left side) to nearly match BFLASH (right 
side) at each data point--the minimum standard deviation in B is auto- 
matically and simultaneously assured. Thus, the exponents on all of 
the terms, other than V , are made negligibly small, and the results 


84 



are trivial, i, e. , the exponents on eq (3-20) are chosen so that eq 
(3-20) is reduced to B % BFLASH. Combining the V and V terms, 

O 

to form the MTWO parameter, avoids this difficulty. Recall that 

0 5 

MTWO oc V (BFLASH) ' , see eq (3-29). Because V is a strong 

o o 

correlating parameter, the computer cannot simply match B and BFLASH 
on either side of eq (3-20); to do so would result in large standard devia- 
tions of B, because the V q component of MTWO would not be properly 
exponentiated. To achieve a minimum standard deviation, as demanded 
in the computer program, the exponent selected for MTWO must be 
compatible with the V and V components of MTWO. 

O 

To substantiate this argument, we refer again to table 5. 3- -note 

that E2 does not vary much when MTWO is substituted for (lines 1 

and 6, 2 and 7, etc. ). This result indicates that V is the dominant 

o 

component in the MTWO parameter,, Also, the mean value of E2, on 

the MTWO term, is ~ 0. 55 (lines 6 to 9)--then the functional dependency 

of B, upon BFLASH, is estimated by combining eq (3-29) and eq (3-27) 

0 28 

to obtain B « (BFLASH) 0 # This functional dependency is well 

removed from the trival relationship, B ^(BFLASH) 

It is remarkable that the variation in the exponent on MTWO is 
so small (lines 6 to 9). Also, the agreement of all of the exponents for 
the hydrofoil and venturi data (specifically lines 8 and 9) is quite good,, 

It appears that the a term could be eliminated in lines 8 and 9. 

The diameter (size) terms in eq (3-12) and eq (3-27) were not 
included in table 5. 3, because our tests were conducted with only one 
hydrofoil and only one venturi size. Moore and Ruggeri [3] obtained an 
exponent value of -0. 1 for the diameter term--E6 = -0 o 1 in eq (3-12)-- 
based on tests using refrigerant 114 in two different venturi sizes. 

Those tests were performed with a venturi identical to the one used in 
our study [20], and with a larger (1. 414:1) geometrically similar venturi. 
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The arithmetic mean value of the developed cavitation parameter, 

K , does not vary appreciably for the venturi data presented in 

c, min 

table 5. 3. This parameter is also relatively constant for the hydrofoil 

data: this is an important result, because constant K . , eq (3-12), 

1 c , mm 

and the isentropic flashing theory [23] are used to predict [4, 5] the 

cavitating performance of a particular piece of equipment. The fact 

that K is different, for different models, limits the current pre- 

c , mm 

dictive techniques [4, 5] to a particular piece of equipment, i. e„ , the 
geometry (shape) of the cavitating equipment must be identical or simi- 
lar. Actually, it was anticipated [20] that K m i n f f° r man Y cavitating 
bodies, would not remain constant--as with the venturi--for all fluids, 
cavity lengths, velocities, temperatures, etc. Then, it was no sur- 
prise that K . for the hydrofoil varied slightly, see table 5, 3. 
c , mm 

Also, K varied more for the hydrofoil than for the venturi; 

c , mm __ 

K was within 7 percent of K . for the venturi, but a similar 

c,mm c, mm 

comparison for the hydrofoil shows 15 percent deviation. An attempt 

was made to mathematically determine the functional dependency of 

K . , upon x, V , and T ; however, the variation in K . was 
c, mm o o c,min 

sufficiently small, and data scatter was sufficiently large--recall that 

instrument error is amplified in this type of parameter- -to discourage 

this effort* It was determined that the hydrofoil K , is independent 

c,mm 

of cavity length, and increases very slightly with increasing tempera- 
ture and velocity, for both hydrogen and nitrogen; the temperature 

effect is stronger than the velocity effect. Because K . varies 

c,mm 

somewhat with flow conditions, for each body, it is convenient to use 

K . for predictive purposes. K . , of course, varies widely 
c,mm ^ c, min 

with body or equipment geometry, as does the pressure coefficient, 

C . Thus, it is quite obvious that prediction of cavitation performance, 

P 

from one piece of equipment to another, will require significant advances 
in the 1 state -of- the art. 1 
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The conventional cavitation parameter for developed cavitation, 

K , also varies with flow conditions for any particular geometry, e. g. , 
see table A- la and Rouse and McNown [28]. 

The results obtained herein indicate that eq (3-27), rather than 
eq (3-12), should be used for predictive calculations [4, 5]. 

The pressure depression in the cavitated region is determined by 
subtracting the measured cavity pressure, in one case, and the satura- 
tion pressure associated with the measured cavity temperature, in the 
other case, from the vapor pressure of the liquid entering the test 
section. In the hydrogen data reported here, the measured cavity pres- 
sure, P , was less than bulkstream vapor pressure by as much as 
1 2 

15. 76 psi (10. 87 N/cm ); these pressure-depressions are obtained by 

subtracting P from P in the tabulated data of appendix A. For the 

1 V 2 

nitrogen data, P exceeded P^ by as much as 10. 08 psi (6. 96 N/cm ). 

Cavity pressure depression increases with increasing cavity 

length, fluid temperature, and velocity for these tests. The data in 

appendix A readily disclose these trends if we permit only one of the 

three parameters (x, V , and T ) to vary at a time and note the value 

o o 

of P - P^ for two different values of the parameter being varied. As 

an example, T and Y are relatively constant in Runs 244 D and 244E, 
o o 

but the cavity length differs by a factor of three; P - P^ for the longer 

cavity (Run 244 E) is about 1. 2 times the P - P, for the shorter cavity 

v 1 

(Run 244 D). Similar comparisons may be made for variations in V 
and T , or portions of the data may be plotted collectively on single 
graphs, etc. , to demonstrate these trends. These functional depend- 
encies are also shown by simply observing the characteristics of the 
experimentally derived exponents in table 5. 3. In reference [23], it 
is shown that the pressure depression increases with increasing T^ and 
B. Referring to line 3 of table 5. 3, we observe that B increases with 
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increasing V , x, and v. Then, for the combined fluid hydrofoil data, 

P - P, must increase with increasing T , V , x, and v. From lines 
v 1 o o 

8 and 9, we observe that P - P, also increases with increasing values 

v 1 

of MTWO. Similar deductions may be drawn for any body-fluid 
combination. 

5. 5 Developed Cavity Shapes 

One of the main objectives of the hydrofoil experiment was to 
obtain cavity volume -thickness data, in an effort to improve the corre- 
lative theory. The hydrofoil-tunnel configuration was designed to pro- 
vide optimum photographs of the developed cavities. Enlarged photo- 
graphs of the cavities, for each experimental data point recorded in 
appendix A, were carefully studied to determine cavity shape, thickness, 
and volume. All of the cavities were elliptically shaped, and the photo- 
graphed cavities were easily fit with a transparent-plastic elliptical- 
template. By recording appropriate data from the template, e. g. , 
major and minor axes dimensions, maximum cavity thickness, and 
angle of projection, it was possible to compute cavity volumes, shapes, 
etc. We found that cavity thickness and volume increased with increasing 
cavity length, and were nearly independent of and T . Because we 
are primarily interested in the shape of cavities near their leading edge, 
we restricted our attention to cavity volumes in the front-half of the 
cavity; in this way, the ill-defined trailing regions of the cavity are 
avoided. The shapes, of all of the cavities, were adequately represented 

by a simple algebraic expression of parabolic form. For hydrogen, we 

03V o 63 

obtained 6=0. 77x ° and for nitrogen 6 = 0 a 44x * , where both 

v v 

6 and x are in millimeters, and x cannot exceed the cavity half-length. 

Interestingly enough, the cavities formed on a wide variety of 
ogival-shaped bodies [28] can also be fit with similar algebraic expres- 

p 

sions of the form 6 = C x , These cavity data, and our hydrofoil data 

v o 


88 



and shape analyses, provided the basis for assuming the existence of 
parabolic shaped cavities in the general analysis of section 3, These 
data also permitted the selection of a mean value for p (~ 0 o 65) in 
section 3. 1 of this report 0 

Because of its application in the pumping machinery field, pressure- 
head has been included in the data tabulated in appendix A. Mathemati- 
cal conversion of pressure to pressure-head merely requires evaluation 
of the liquid density at the point of measurement; however, selection of 
the appropriate liquid density can be a bit perplexing. Figures 5. 5 to 
5. 13 indicate that the measured pressures and temperatures, within the 
cavities, are not in perfect agreement. Also, due to the thermal expan- 
sivity of liquid hydrogen, the bulkstream temperature does not remain 
perfectly constant as the liquid flows over the hydrofoil. The following 
methods were used to calculate pressure-head from the cavity measure- 
ments: (1) Head (h ) was calculated from measured cavity pressure by 

n 

using the saturation density at the measured pressure, (2) Head (h ) 

n, T 

was calculated from measured cavity temperature by using the saturation 
density at the measured temperature. Both values of head are given in 
the tabulated data in appendix A. 

6. CONCLUDING REMARKS 

Desinent cavity data, for a 0. 5 -caliber hydrofoil, were acquired 
for vaporous hydrogen and nitrogen cavities; the results are given in 
appendix A and on figures 5, 1 to 5. 4. Correlation of the desinent data 
is treated in appendix C. The desinent data tend toward a single-valued 
K , at the maximum velocities, irrespective of temperature. The 
hydrogen data are more temperature -dependent than the nitrogen data. 
Also, these cryogenic liquids require less subcooling--relative to 
higher boiling-point liquids--for desinent cavitation to occur. 
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Pressure and temperature profiles were measured within fully 

developed, vaporous hydrogen and nitrogen cavities; these results, for 

the hydrofoil, are given in appendix A and on figures 5„ 5 to 5, 13. Within 

data accuracy, these pressure and temperature depressions were in 

stable thermodynamic equilibrium. These data were correlated using 

a previously described [20] technique, and the extended theory developed 

herein. Using the conventional correlating technique [20], it was found 

that u , V , x, and v were valuable correlating parameters for combined 
o 

fluids. Using the new MTWO parameter, only o', MTWO, and x were 
of value--it seems quite likely that o' may be negligible for combined 
fluid correlations. Thus, only MTWO and x may be required for future 
correlations. If the MTWO correlation is not used, the results may be 
degraded by ~ 50 percent (as based on standard deviation in B). Because 
MTWO is such an influential parameter, its use is highly recommended 
in future work, for both correlative and predictive purposes. 

Correlative expressions, predicted from convective heat transfer 
theory, and evaluated from experimental data, are in reasonable agree- 
ment; however, more analysis is required to explain the intricacies of 
the mass transfer phenomena- -an important physical effect, as evidenced 
by the influence of the MTWO parameter. 

K was found to vary by about 15 percent, about a mean value, 

c,min 

K . , for the hydrofoil. This was about double the deviation experienced 
c, mm 

with the venturi tests [20]. Variations, in K . , of 1. 65:1 were found 

L J c , mm 

in tests on zero-caliber ogives [24] in water. It is quite apparent that 
K . will vary with equipment geometry, fluid, velocities, tempera- 
tures, etc. Then, the current predictive technique [4, 5] , that relies 

on constant K . (or K . ), must be used with appropriate caution, 
c, mm c, mm 
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The cavity-shape data, acquired during this study, indicates that 

the cavities can be described by a simple expression of the form 

6 = C ; this expression is valid only in the frontal regions of the 

vo 

cavity* 

7. NOMENCLATURE 


a 


l 


a 

v 


B 


BFLASH 


C 

n 


C 

P 

P 



C 

P> v 



acoustic velocity of liquid, evaluated at the minimum 
cavity pressure, 

acoustic velocity of vapor, evaluated at the minimum 
cavity pressure, 

ratio of vapor to liquid volume associated with the 
sustenance of a fixed cavity in a liquid 

B derived from isentropic flashing theory (Ref* [23]) 

BFLASH 

(n = 0, 1 , 2, 3 * * * etc* ): constants or numerical 
coefficients in various algebraic expressions 

2 

pressure coefficient [ = (h - h )/(V /2g )] 

L ~ x o o to c 

v Z 

minimum pressure coefficient [ = (h - /2g^)] 

specific heat of liquid at constant pressure, evaluated 

at the tunnel inlet (P and T ) 

o o 

specific heat of vapor at constant pressure, evaluated 
at the minimum cavity pressure, P^ 

specific heat of liquid at constant volume, evaluated 
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C = specific heat of vapor at constant volume, evaluated 

v, v 

at P 

D = characteristic dimension of cavitating body, e. g. 

diameter of an axisymmetric body, or thickness (t) 
of a hydrofoil 

g = conversion factor in Newton's law of motion (gravita- 

c 

tional acceleration) 


G 

c 


h 

n 



h 

o 


h 

v 


choking two-phase liquid-vapor mass -flow flux, see 
eq (3-21) 

(n = 1, 2, 3, 4, or 5): head corresponding to cavity 
pressure, measured at a particular instrument port 
on the hydrofoil 

(n = 1, 2, 3, 4, or 5): head corresponding to the 
saturation pres sure at the cavity temperature, 
measured at a particular instrument port on the 
hydrofoil 

tunnel inlet head corresponding to absolute inlet 
pressure 

head corresponding to saturation or vapor pressure 
at the tunnel inlet temperature 


= convective heat transfer coefficient at axial distance x 

= head corresponding to absolute pressure, measured 

on the hydrofoil at distance x, downstream of the mini- 
mum pressure point- -for non- cavitating flow 


92 



V 

h 


k 


K 

c, min 
K 

c, mm 


K. 

IV 


K 

v 


M 


2 cp 


MT WO 


Nu 

x 

Pe 

P 

n 



P 

o 


= head corresponding to the minimum absolute pres- 
sure on the leading edge of the hydrofoil, computed 

v 

from expression for C 

P 

= thermal conductivity of liquid, evaluated at tunnel 
inlet 

= developed cavitation parameter , based on minimum 

2 

measured cavity pressure [ = (P - P,)/(p V / 2g )] 

o 1 o o °c 

= arithmetic mean value of K . for a complete set 

c, min 

of data points for a particular hydrodynamic body- 
fluid combination 

= cavitation parameter, K , evaluated at incipient 
(desinent) conditions [ = (P q - P ) /( /2 g^)] 

2 

= developed cavitation parameter [ 5 (P^ - P )/(P q V o / 2 g^)] 

= two-phase Mach number [ = V /V_ 1 

c L ~ o 2 cp 

= liquid phase velocity ratio [ = V /V ] 

O i/ 

= length Nusselt number [ = ft^x/k] 

= Peclet number f= V x/ a] 

L o 

= (n = 1, 2, 3, 4, or 5): absolute cavity p re s sure , 

measured at a particular station or instrument 
port on the hydrofoil 

= (n = 1, 2, 3, 4, or 5): saturation pressure cor- 

responding to the measured cavity temperature at a 
particular station or instrument port on the hydrofoil 

= tunnel absolute inlet pressure 
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Pr 

P 

v 

R 

Re 

t 

T 

n 


T 

o 



Prandtl number [ = C g,/k] 

P> ^ 

saturation or vapor pressure at tunnel inlet 
temperature 

radius of axi symmetric body, see figure 3. lb 
Reynolds number [= pV^x/p,] 
maximum thickness of hydrofoil 

(n = 1, Z, 3, 4, or 5): measured cavity temperature 
at a particular station or instrument port on the 
hydrofoil 

bulkstream temperature in degrees Rankine (Kelvin), 
of liquid entering the tunnel 

characteristic liquid velocity component, normal to 
cavity wall 


V 

o 



velocity of test liquid at inlet to tunnel 

velocity of two-phase liquid-vapor mixture across 
cavity wall 



-v 

v 


We 

x 


X 


volume of liquid associated with the sustenance of a 
developed cavity 

volume of vapor associated with the sustenance of a 
developed cavity 

Z 

Weber number [ 5 pV x/a] 

axial distance measured from minimum pressure 
point on hydrofoil, used to measure cavity length, 
and as an axial coordinate in analyses 
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X 


z 

Greek 

a 



a 


fluid quality, mass fraction of vapor per unit mass 
of fluid 

evaluated by eq (3-23) 


thermal diffusivity of liquid, evaluated at tunnel inlet 

characteristic thickness of the thermal boundary 
layer in the liquid 

thickness of the developed vaporous cavity 

absolute viscosity of liquid, evaluated at tunnel inlet 

kinematic viscosity of liquid, evaluated at tunnel 
inlet [ E mV P] 

density of liquid, evaluated at tunnel inlet [ = p^] 

density of liquid, evaluated at minimum cavity 
pressure, 

density of vapor, evaluated at minimum cavity 
pressure, P^ 

density of two-phase liquid-vapor mixture at the liquid- 
vapor interface, evaluated at P^ 

surface tension of liquid in contact with its vapor, 
evaluated at tunnel inlet 


Subscripts 

l = denotes liquid phase 

o = denotes tunnel inlet location 
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ref 


v 

2cp 

Superscripts 

El 


E2 


E3 


E4 


E5 


E6 


E7 

ml 


reference run (data point), or test conditions, to 
which a computation is being referenced when 
attempting to correlate cavitation performace via 
eq (3-12) or eq (3-27) 

denotes vapor phase 

denotes two-phase liquid-vapor 


exponent on thermal diffusivity ratio in eq (3-12) 
and eq (3-27) 

exponent on tunnel inlet velocity ratio in 

eq (3-12) and also used as an exponent on the MTWO 

ratio in eq (3-27) 

exponent on cavity length ratio in eq (3-12) and 
eq (3-27) 

exponent on kinematic viscosity ratio in eq (3-12) 
and eq (3-27) 

exponent on surface tension ratio in eq (3-12) and 
eq (3-27) 

exponent on (characteristic dimension) cavitating 
body diameter (or thickness) ratio in eq (3-12) 
and eq (3-27) 


exponent on ratio in eq (3-20) 

exponent on Reynolds number in Nusselt relationship, 
see eq (3-1) 
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exponent on Reynolds number in thermal boundary 
layer expressions, see eq (3-2) 

exponent on Prandtl number in Nusselt relationship, 
see eq (3-1 ) 

exponent on Prandtl number in thermal boundary layer 
expressions, see eq (3-2) 

exponent in algebraic expression for cavity shape 
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APPENDIX A: Experimental cavitation data— nitrogen and hydrogen— for hydrofoil. 
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Table A- lb. (cont'd) 
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Table A-2a. Experimental cavitation data for 0.312-inch hydrofoil using liquid hydrogen (English Units). 
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Table A-2b. (cont’d) 
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Table A-2b. (cont’d) 


I- 

00 

o 

o 

4 

in 

r- 

O' 

CO 

CO rH 

O' rH 

CM 



On O 

O' 4 


s0 

r- 

nO in 

o 

n- 


« 

• 

• 

• 

• 

• 

• 

• 

• • 

• • 

• 

• 

• 

• • 

• • 

• 

• 

• 

• • 

• 

• 

i n X 

4 

in 

CM 

m 

O' 

rH 

On 

vO 

M3 CM 

r- cm 

m 

o 

o 

CO CO 

4- 4 

On 

CO 

M> 

CO o 

m> 

o 


H 

eg 

co 

CO 

00 

o 

CM 

n 

in n 

rH in 

o 

in 

M3 

O CM 

m O' 

in 

M) 

rH 

M3 4 

ao 

r- 

X 

rH 

t-» 

rH 

rH 


<— 1 

rH 

rH 

rH rH 

rH rH 

rH 

rH 

rH 

CM <M 

rH 

rH 

fH 

rH 

rH CM 

rH 

rH 


1- 

CM 

fH 

o 

e- 

St 

O 

rH 

CM 

O' v0 

CO O 

r- 

ON 

m3 

O CO 

CM 4- 

in 

4 

co 

CO M> 

O 

CO 

* 

a 

• 

• 

• 

• 

• 

• 

• 

• • 

♦ • 

• 

• 

• 

» • 

• ♦ 

• 

• 

• 

* • 

• 

♦ 

4 X 

CM 

O' 

CO 

CO 

rH 

4 

4 

CO 

CM O 

r- co 

CO 

CO 

m 

O' rH 

CM sO 

CO 

n 

m 

rH 4 

O' 

O' 


rH 

CM 

CO 

CM 

O' 

O' 

CO 

rH 

4* CM 

On in 

On 

in 

4 

N CO 

CM r- 

m 

in 

O' 

CM CO 

r- 

CO 

X 

rH 

rH 

rH 

rH 



rH 

rH 

*H rH 

fH 


rH 

rH 

•H CM 

rH rH 

rH 

*H 


rH CM 

rH 

rH 


t- 

O' r- CO CM 

o cm r- 

00 

M) 1^ 

CO M3 

f— 4 oo 

4 M3 

M3 rH 

rH O rH 

M3 co 


in 

CO X 

CM *0 r- CM 

r*- o co 

co 

r- 4- 

CM h- 

co co 

cnj in 

nO 4 

in in rH 

On O' 

•H 

00 


H H(\| H 

r- oo h 

o 

CO O' 

Q0 CO 

CO 4 CM 

4 rH 

O' CO 

4 r0 O' 

O 00 

r- 

o 

X 

rH rH rH rH 

rH 

rH 

iH 

rH 

rH fH 

H CM 

rH 

fH rH 

rH rH 

iH 

fH 


1- 

\0 o o m 

ifUTlN 

h- 

CM v0 

I s - 00 

4 H co 

O' CM 

r- 4 

4 M3 O 

co r- 

00 

in 

CM X 

vO 4 rH co 

00 O rH 

NO 

4 in 

h- fH 

rH 4 O' 

CO CO 

no in 

n t in 

00 CM 

00 

r- 


O' O CM O 

N 00 O 

■H 

CM O 

00 CM 

O' co <h 

CO O' 

o co 

CM CM O' 

cm r- 

4 

fH 

X 

fH rH rH 

rH 

rH 

rH iH 

rH 

rH rH 

H i — 1 

rH fH 

rH rH 

rH rH 

fH 

rH 


*- 

co co 4 O 

CO O' CO 

CO 

r- r- 

CO CM 

O CNJ 00 

o r- 

CO 00 

03 co r- 

O' CO 

O' 

O' 

rH X 

m m r- r- 

O' O' 4 

O 

CO O' 

4 cm 

r- O H 

S rH 

00 CO 

4 CO rH 

r- r- 

O' 

O' 


O' O' O O' 

f- If On 

rH 

•— 1 O' 

00 fH 

CO CM rH 

CM sO 

O' CM 

rH rH O' 

rH m 

CM 

o 

X 

rH 


rH 

rH 

rH 

rH tH 

rH rH 

rH 

H rH 

fH rH 

rH 

rH 



rH CM 

00 

CO 

M3 

M> 

NO 

CO 

rg rH 

CO vO 

Ch 

CO 

O' 

00 00 

O in 

NO 

O' 

CO 

CO oo 

CM 

o 

in 

• • 

• 

* 

• 

• 

• 

• 

* • 

• • 

• 

• 

• 

• • 

• • 

• 

• 

• 

• • 

• 

• 

X 

NO O 

in 

oo 

rH 

CM 

in 

fH 

CO O' 

rH 4 

4 

CM 

o 

If' CO 

00 sO 

CO 

in 

CO 

cm in 

CO 

O' 

X 

CO NO 

M3 

4 

O 

CM 

e- 

nO 

M3 NO 

4 m 

O' 

rH 

O' 

O 4 

4 O' 

4 

co 

o 

m o 

4 

oo 


fH rH 

fH 

rH 

rH 

rH 

iH 

rH 

CO CM 

rH CM 


CO 

CM 

CO CO 

CM CM 

CM 

CM 

H 

•H 4 

CO 

CM 



CM O' O rH 

4 O' nO 

00 

in rH 

co m 

CO 00 4 

o r- 

M3 M3 

in r- in 

CO O' 

rH 

fH 

X 

00 4 O O 

O' 4 

00 

fH O' 

M3 CM 

MD CO fH 

00 o 

4 00 

O in rH 

O c\j 

o 

CO 

X 

o in so co 

oo O' m 

00 

4 rH 

GO 4 

I s - O O 

CM CO 

rH CO 

CO 00 CO 

rH M3 

CO 

4 


rH rH rH i-H 

rH 


CO rH 

CM 

CO CM 

CM CO 

fH rH 

CM rH 

rH CO 

CO 

i-H 


co 

o co oo r- 

00 fH co 

00 

4 O' 

00 rH 

O CM M> 

co m 

h m 

4 0 0 

CM O' 

O' 

in 

X 

rH o in oo 

CD 00 rH 

r- 

O O' 

r- cm 

nO CM 4 

M3 4 

in o 

O oo oo 

O' O 

in 

o 

X 

O' 4 in O' 

00 00 rH 

r- 

co r- 

r- cm 

e- cm O' 

CO 00 

O' rH 

co O' h- 

O M3 

CO 

o 


rH i-H 

rH 


r-H 

iH 

CM 

rH CM 

r— 1 

fH 

rH rH 

CM 

rH 



4 M3 f'- 

rH 

CD CM O 

oo 

fH O' 

oo m 

CM 

O' 

M3 

4 00 

-4- CO 

O' 00 o 

r- r- 

CO 

r- 

X 

O' CO 00 

CM 

r- O' r- 

in 

CM M3 

rH sO 

O' 

r 

O 

rH rH 

4 co 

CM r 03 

r- m 

00 

CO 

X 

00 O rH 
rH rH 

O' 

CO CO o 

00 

O CD 

rH 

CO O' 

h- 

o 

iH 

O' 

CO s0 
rH fH 

O' o 

rH 

O' 00 N 

O 4 

rH rH 

CM 

r—i 

o 

rH 



O' O CO so 

oo oo 

o 

CO 

M3 O 

OO CO 

CO 

CO fH 

4 CM 

in o 

CO CO o 

r- O' 

O 

00 

X 

O' r- o o 

r- oo 

r- 

n- 

m rH 

CO oo 

o 

4 in 

rH O' 

O' rH 

oo r- oo 

O' O' 

O 

H 

X 

00 O' O O' 

co 00 

O' 

00 

O O' 

00 O' 

00 

O O' 

CO 4 

O' rH 

co co r- 

O 4 

CM 

i-H 


rH 




rH 



fH 

rH rH 

<H 


rH rH 

rH 

rH 



DUO 

UJ 

DUO 

< 

CD U 

co u 

CD U C 

CO U 

CO u 

DUO 

OJ u 

< 

03 

x • 

4 4 4 

4 

in m in 

M3 

r- r- 

00 00 

O' O' O' 

fH f — 1 

CM CM 

4 4 4 

in m 

M3 

I s - 

3 O 

4 4 4 

4 

4 4 4 

4 

4 4 

4 4 

4 4 4 

in in 

m m 

m in in 

m m 

m 

m 

CL Z 

CM CM CM 

CM 

CM CM CM 

CM 

CM CM 

CM CM 

CM CM CM 

CM CM 

CM CM 

CM CM CM 

CM CM 

CM 

CM 




143 



b- 

M> O' 

CM rH 00 

cn <r r- 

<M \0 nj 

O CO 

f- cm 

r- 0 n 

m 

4 O' 

O 

O' 

tn X 

in 4 

h cn cn 

O' cm 4 

0 n* O 

in co 

r- M3 

CO 03 00 

0 

CM CO 

CM 

CO 


vO O' 

co m> r- 

co 4 4 

cm co in 

cm r- 

0 n 

CM rH CO 

CM 

03 CO 

CO 

O' 

X 

rH rH 

HHrt 

rH rH rH 

«H H rH 

rH rH 

CM CM 

CM CM cm 

CM 

CM CM 

CM 

CM 


h- 

O 

in CM I s - 

CO co in 

CO CO CM 

O' O' 

O' CM 

rH CO 

tn 

M3 O 

r- 

rH 

4 2 

in 00 

O O 

0 cn H 

-4" co co 

cm n 

CO sO 

00 00 CO 

an 

CO CM 

0 

rH 


co 0 

rH in I s - 

(O CO 4 

rH in 

rH r- 

in in 

CM CM CO 

0 

r- r- 

rH 

CM 

X 

rH rH 

H H H 

rH rH rH 

rl H rl 

rH rH 

rH CM 

CM CM CM 

CM 

CM CM 

CM 

CM 



in cm 

O' r-~ cm 

CO CO 0 

CM 03 CO 

cm in 

M) O' 

co M3 4 - 

4- 

in cn 

O' 

'O 

CO 2 : 

co CO 

CO M) CM 

CO M3 CM 

O H CO 

O' CO 

CM 00 

O' rH O' 

CM 

vO co 

O' 

4 


O 00 

to h r 

O rH 4- 

O' CM 4 

00 co 

CM O 

4- O' CM 

CO 

CO O' 

co 

M3 

I 

rH rH 

rH rH 

rH rH r-t 

rH *— 1 

rH 

rH CM 

H HIM 

rH 

CM rH 

H 

|H 


b- 

CO M3 

4- CO 00 

h in 0 

O CM CM 

4* 00 

IH O' 

00 h 4 

CO 

co 00 

M3 

CM 

CM X 

rH CM 

CM 4 M3 

co co CM 

r- in 0 

rH I s - 

co 4 

I s - oo m 

m 

00 <0 

M> 

rH 


•“» M3 

O' rH in 

O 0 co 

0004 

O' CM 

CO 00 

in M3 0 

4 

O O' 

M3 

00 

X 

rH rH 

iH rH 

rH rH rH 

rH rH 

rH 

rH rH 

rH rH CM 

rH 

CM rH 

rH 

rH 


(- 

CM 00 

4 00 O' 

4 4 cm 

O' I s - n 

4 m 

00 rH 

O' M3 n- 

4 

CM O' 

00 

rH 

rH X 

in co 

rH r- CM 

rH CO 4 

rH 4 CM 

co r- 

in in 

00 n O' 

in 

rTn 

00 

M3 


0 4 

O' 0 (O 

O' O' 0 

OO O' rH 

O' rH 

CM M3 

<0 4 M3 

CM 

00 I s - 

4 

M3 

X 

rH <H 

rH rH 

rH 

rH 

rH 

rH rH 

rl rl rl 

rH 

rH rH 

rH 

rH 



*h ur\ in m tn n- m o 4 o co o <j- co cm rH o> on in oo 

in •• ••• • * • ••• •• •• ••• • •• • • 


X 

rH 

cm in co 

O CO M 3 

in n 00 

O' 4 

O' M 3 

CM rH O' 

4 

rH 4 

0 

M3 

X 

0 0 

co 4 n 

M 3 r- 00 

M3 00 O' 

CM CO 

rH 4- 

M 3 OO O' 

CM 

n cm 

0 

I s - 


co 4 

H CM CM 

rH rH rH 

rH rH rH 

rH CM 

CM CO 

CM CM CM 

CO 

co <0 

CO 

CO 


~Q 




00 O' 

00 CM CM 

CD CO O' 

.7 

,7 

0 

M3 in 

CM O 

r- 4 in 

M3 

O O 

r- 

r- 

* x 

r- m> 

O M> CO 

in 00 0 

4 co 

rH 

4 M3 

00 in 

O co O' 

00 

rH 4 

CO 

M3 

X 

in ao 

rH in M3 

cm tn 00 

CM I s - 

O' 

0 n 

CO CM 

CM M3 00 

I s - 

cn r- 

n* 

rH 


rH CO 

rH rH CM 

r-i r—i r^-i 

ft rH 

rH 

rH CM 

rH CO 

(M CM CM 

rH 

CO CM 

CM 

CM 


3 

M 

C0 «H 

O' 00 r- 

M3 4 O 

O O' O' 

rH 00 

co r- 

CM 

CM 4 

00 

4 O' 

M3 

r- 

4 M3 

CM CO CO 

M3 r- in 

00 4 O 

CM rH 

co r- 

O' 

in cm 

CO 

'0 0 

rH 

<0 

X 

O CM 

O' O CO 

O' O M> 

O' CM 00 

O co 

cm h* 

4 

00 in 

4 

cm r- 

I s - 

sO 


rH CM 

rH CM 

rH rH 

rH rH 

rH rH 

rH rH 

rH 

rH CM 

iH 

CM rH 

rH 

r—i 



4 in 

O' CM H 

in M3 00 

4 <0 00 

«H O 

00 n 

n 4 0 

4 

I s - 4 

CM 

CO 


00 co 

CO O »H 

rH 4 rH 

4 CO I s - 

CM O 

I s - M3 

4 4 00 

in 

O CM 

CO 

in 

X 

O CO 

O' O CO 

O' O' rH 

O' O CM 

O CM 

cm in 

4 n I s - 

4 

CO M3 

n 

M3 


rH rH 

rH rH 

rH 

rH rH 

«H rH 

rH rH 

rH rH rH 

rH 

t*H pH 

1 — 1 

*H 



m O' 

O CM CO 

4 rH in 

O' r- cm 

M3 O 

O' O' 

0 00 n 

0 

4 4 

in 

4 

^ X 

rH rH 

M3 O' O 

00 rH 00 

rH M3 M3 

rH O 

rH r- 

M3 CM CM 

I s - 

n ih 

e 

I s - 

X 

rH CO 

O' O' rH 

00 O' O' 

O' O' 0 

O CM 

co in 

4 m M3 

4 

I s - M3 

n 

M3 


rH rH 

rH 


rH 

rH rH 

rH rH 

«H rH rH 

rH 

rH rH 

rH 

<H 



ac 0 

an u a 

cd u a 

QJ U 6 

CD U 

00 u 

CO u O 


CD ^ 

CD 

< 

Z • 

co 00 

000 

<H rH rH 

CM CM CM 

CO CO 

4 4 

in in n 

O' 

O O 

CO 

rH 

3 O 

in in 

M3 M3 M3 

M3 M3 M3 

M3 M3 M3 

M3 M3 

M3 M3 

M3 M3 M3 

M3 

h- I s - 

r- 

00 

or z 

CM CM 

CM CM CM 

CM CM <M 

CM CM CM 

CM CM 

CM CM 

CM CM CM 

CM 

Cvj CM 

(M 

CM 


144 



Appendix B: Thermocouple fabrication procedure 

After the thermocouple sensor ports are drilled into the assembled 
hydrofoil- -see figures 4. 2 and 4. 5- -the thermocouples are fabricated 
and installed as follows: 

I. Deburr thermocouple sensor ports 

A. Thread a 0. 010-inch (0. 0254 cm) piano wire through the 
sensor port and through the small stainless steel conduit 
tube. 

B. Apply a medium- coarse lapping compound to the wire at the 
sensor port. 

C. Slide the piano wire back-and-for th inside the conduit ~ 200 
times to remove any burrs, at the sensor port, that may 
have resulted from the drilling operation. 

D. Repeat this deburring operation for each thermocouple port. 
When the last port has been deburred, remove the piano wire 
and flush all conduits with solvent. 

E. Rod dead end of each conduit with piano wire and flush with 
solvent. Repeat this operation, as required, until thermo- 
couple conduits are clean. Dry conduits by blowing com- 
pressed air through them. 

II. Fabricate thermocouple junctions 

g 

A, Cut a length of Chromel wire twice the length required to 
extend through the conduit. Skin the thermocouple wire 
insulation at one end and in the center; tin the skinned wire 
using stainless steel soldering flux and solid core (50/50) 
solder. Wash off flux. 


Identification of a manufacturer and a manufacturer’s product has been 
necessary to make the results of this work meaningful and in no way 
implies a recommendation or endorsement by the National Bureau of 
Standards, or by the National Aeronautics and Space Administration. 
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B. Cut a length of Au (0. 07 atomic percent Fe) wire one-half 
as long as the Chromel wire, and skin and tin each end. 
Use voltage-controlled soldering iron and do not get Au 
wire too hot- -use just enough heat to melt solder. 

C. Twist the Chromel-Au wires together to form a junction, 
at the mid-length of the Chromel wire, and solder using 
resin-core radio solder. Clean junction with solvent and 
thoroughly dry junction. 

D. Reinforce the enamel insulation on the Chromel-Au wires 
with polyimide enamel--two coats. Wipe the enamel on 
with a plastic glove. Do not insulate the junction (stop 
application of enamel at about 0. 6 cm from the junction). 
Keep the Chromel and Au wires separated until the enamel 
dries. 

Installation of junctions 

A. Thread the insulated end of the extra-long Chromel wire 
into the tube -conduit from the sting end of the hydrofoil 
assembly. The extra length of Chromel wire is used as 
a fish-wire to pull the thermocouple junction through the 
conduit. The thermocouple junction is pulled through the 
sensor port so that all of the uninsulated Chromel-Au 
extension wire, near the junction, is exposed. 

B. Apply polyimide enamel to the uninsulated wires and 
junctions --two coats --allow the enamel to thoroughly dry. 

C. Pull the junctions back into the sensor ports (final desired 
position is flush with hydrofoil surface) and check for 
electrical shorts. 

D. Attach a transparent vacuum hose to the aft end of the 
hollow cylindrical sting shaft. 


9 

E. Mix epoxy: Emerson and Cuming, Inc. , Stycast 2850 GT 
with room temperature cure activator. 

F. Apply epoxy to the sensor ports (being careful not to displace 
the thermocouple junctions); then use the vacuum assist and 
carefully manipulate the thermocouple wires until epoxy can 
no longer be sucked into the conduits through the sensor ports. 

G. Remove vacuum hose and excess epoxy from each sensor 
port; leave the epoxy at each of these ports slightly mounded. 

H. Make final placement of the thermocouple junction at each 
sensor port. Junctions should project slightly from the 
mounds of epoxy on top of each sensor port, and should be 
visibly centered in the sensor ports; i. e. , the thermocouple 
junctions should not be in contact with the hydrofoil sensor 
ports. These soldered junctions should be sufficiently long 
that they will not be destroyed, later on, when the epoxy 
mounds, and encapsulated junctions, are filed and sanded 
flush with the hydrofoil. 

I. After the epoxy is thoroughly cured, the thermocouple instal- 
lation is again checked for electrical shorts. Then tape-mask 
the hydrofoil, exposing only the epoxy bumps. Use flat Swiss- 
pattern files and file the bumps and encapsulated junctions 
until they are nearly flush with the hydrofoil. Then use wet 
fine -grit sandpaper and sand the epoxied junctions flush with 
the hydrofoil. Again check for electrical shorts. These 
sealed and flush mounted junctions are shown in figure 4. 1. 


Identification of a manufacturer and a manufacturer's product has 
been necessary to make the results of this work meaningful and in no way 
implies a recommendation or endorsement by the National Bureau of 
Standards, or by the National Aeronautics and Space Administration. 
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J. Wire motion, and consequently chafing, is minimized by 

pouring an epoxy ping in the conduit support assembly at the 
aft end of the hollow sting shaft- -this minimizes the possi- 
bility of electrical shorts developing. The tinned Chromel- 
Au wires extend to the rear of this epoxy plug--see figure 4. 2. 

IV. Thermocouple extension wires 

A. Thermocouple extension wires are then attached to the 
tinned junction leads, at the rear end of the sting assembly, 
using radio solder. These connections are then enamel 
insulated, and the extension wires are routed through one 
of the hollow elliptical sting struts, and soldered to a gold- 
pin electrical connector--see figures 4. 1 and 4. 2. 

B. About 75 cm of each extension wire is bundled and carefully 
coiled and compactly packed with a stiff vacuum grease, to 
form an ogival shaped projection, on the rear end of the 
sting shaft. The hollow ogival sting cap- -see figure 4. 2-- is 
then screwed onto the sting shaft, completely protecting the 
thermocouple wires. The vacuum grease packed wire 
assembly provides some thermal tempering for the thermo- 
couple extension wires. 

C. A low temperature vacuum seal is obtained at the gold-pin 
electrical connector by casting another epoxy plug. 

D. Constant temperature at this electrical connection was 
ensured by 1) using a low resistance (brass) thermal path 
between the connector and the sting flange assembly and 2) 
carefully insulating this connector with multilayer radiation 
insulation. 
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Appendix C: Computer program for correlating 

desinent cavitation data 

This appendix outlines the details of the computer program used 
to correlate the desinent and incipient cavitation data acquired from 
this study. Because data could not be obtained at constant temperature, 
a mathematical technique was developed, so that isotherm data could be 
calculated. While the fluid temperature could not be precisely controlled, 
the data were obtained in discrete groups that were nearly isothermal; 
thus, the technique for computing isotherm data need not be too 
sophisticated. 

The technique chosen to correlate these desinent data consists of 
fitting an equation to a three-dimensional surface, i. e. , P , V , and T^ 
(at desinence) are the three coordinates, and an equation is derived to 
fit all of the experimental data. This equation may then be used to pre- 
dict desinent flow conditions (P and V ), along any isotherm (T = 

o o o 

constant), within the boundaries of the experimental data. Also, the 

desinent flow conditions, at each individual experimental data point, may 

be computed for comparison with each actual experimental data point. 

Actually, this is the technique used to assure that the equation adequately 

fits the experimental data, i. e. , the equation is derived from a least- 

squares solution technique that minimizes the differences in calculated 

P , and experimental P , at each experimental T and V . When these 
o o o o 

differences in P^ are each less than the minimum possible experimental 
error in P^, the equation is considered adequate. 

The computer program used to derive this correlative equation 
uses a simplified Bjorck subroutine. This Bjorck technique is pro- 
grammed to provide a conventional least-squares solution to data-fitting 
problems. This program always locates the absolute minima in the 
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selected ordinate--in this case, standard deviation in P at desinence, 

o 

The Bjorck program consumes little machine time and is very reliable. 
More details concerning the application of the Bjorck program, and its 
niceties, to data-fitting problems, are given in reference [20], 

To enter the Bjorck program, it is necessary to assume an appro- 
priate form for the surface -fitting equation. Previous experience with 

P - V - T data, and preliminary hand calculations, indicated that a 
o o o 

simple polynomial expression would suffice. The general form of this 
expression is 


P = A T + A V + A T V +A.T 2 + A V 2 
o lo 2o 3oo 4o 5o 


4- A,(T V) + + A , 

6 o o n 


(C-l) 


where the A to A are numerical coefficients. It was soon learned 
1 n 

that the first four or five terms, of eq (C-l), provided an adequate fit 
of the experimental data and the higher order terms were eliminated 
from the program. The program listing on page C-6 shows that only 
eight terms were retained in the program; the minimum number of 
these terms were used to obtain a satisfactory fit of each batch of 
experimental data. 

As expected, the order in which the terms of eq (C-l) are intro- 
duced into the Bjorck program influences the number of terms required 
to obtain a satisfactory fit of the experimental data. Thus, various 
combinations of terms--where each combination contained not more 
than eight terms--were tested to obtain the best combination for each 
batch of data. The Bjorck solution will best fit the data, with the mini- 
mum number of terms, when these terms are the most influential of 
the available terms, and appear first in the mathematical arrangement 
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of eq (C-l). The best combination and number of terms is readily 

indicated, by the minimum standard deviation in P^, i. e. , the Bjorck 

program computes this standard deviation fori, 2, 3, 4, 5, 6, 7, and 

8 terms, for each combination of terms. Once the best arrangement 

and number of terms are selected, the differences in P (calculated 

o 

minus experimental values) are inspected to see that the data have been 

fit, within experimental uncertainties. If so, this correlative equation 

is used to predict P , at constant T , for various values of V , within 

o o o 

the boundaries of the experimental data. This equation should not be 
used to predict desinent cavitation conditions, outside of the confines 
of the experimental data. As in all surface-fitting cases, it is dan- 
gerous to attempt to extrapolate an experimentally defined surface. 

If the correlative equation does not fit the experimental data, 
within experimental uncertainties, it is necessary to try a different 
combination of terms in eq (C-l), or assume a different general form 
for the correlative expression. The polynomial expression in eq (C-l) 
was sufficient for all data acquired in this study. 

For desinent nitrogen cavities on th§ -hydrofoil, we obtained 

P = -1. 34832 T + 0.81918V - 0.00911 TV 

o o o o o 

+ 0. 01082 T 2 + 0. 01886 V 2 . (C-2) 

o o 

Although incipient cavity data are included in the tabulations of appendix 
A, only desinent data were used in deriving eq (C-2). The incipient 
data were not included because they exhibited a slight hysteresis effect, 
relative to the predominantly desinent data. 
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For desinent (and incipient) hydrogen cavities on the hydrofoil, 
we obtained 

P = -3. 44809 T + 0. 57330 V - 0. 00488 T V 
o o o o o 

+ 0. 08066 T 2 . (C-3) 

o 

The incipient cavity data, for hydrogen, showed negligible hysteresis 

and were included in deriving eq (C-3). The data shown on figures 5. 3 

and 5. 4, and in table 5. 2, were computed from eq (C-2); likewise eq 

(C-3) was used to produce figures 5. 1 and 5. 2, and table 5. 1. In eq (C-2) 

and eq (C-3) the units of P , V , and T are psia, ft/s, and degrees 

o o o 

Rankine , re spe ctively. 

The Bjorck desinent-data fitting program functions as follows: 

(1) A particular combination of polynomial terms is chosen- - 

statement 19; all input data, label (run no. ), T , V , and 

o o 

P are read and stored in computer memory- -statement 10. 
o 

(2) Coefficients for the Bjorck subroutine are computed-- 
statements 1 1 to 14. 

(3) The chosen algebraic terms are used to fit the experimental 
data by calling the Bjorck subroutine --statement 22 + 2. 

(4) Standard deviation and residuals are computed and printed- - 
statements 22 + 3 to 24 + 1. 

(5) The number of algebraic terms to be used in the correla- 
tive expression are predetermined, along with ranges and 

intervals for T and V ; the coefficients for the chosen 
o o 

correlative expression are printed- - statements 30 to 33. 

(6) The differences in calculated and experimental P , at each 
data point, are computed and printed- -statements 44 to 47. 
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(7) 


Desinent pressures (P ) corresponding to the chosen 

isotherms (T ), at the chosen velocities (V ) , are then 
o o 

computed and printed- - statements 47 + 1 to 57 + 9. 

(8) Values of K. , at each of the coordinate positions computed 

iv 

in step (7), are calculated and returned upon demand-- 

statement 30 and statements 9 to 90. Metric (SI units) 

values of P , V , and T , at each of the calculated 
o o o 

coordinate positions, are also computed and returned- - 
statements 61 to 62 + 4. 

All pertinent data is printed out as it is generated. Different 
combinations of polynomial terms, as listed in statements 12 to 14, 
may be tested, on each batch of experimental data, by cycling the 
program with appropriately punched data process cards--statements 
90 and 19. Use of any linear expression, different from eq (C-l), is 
made possible by simply replacing the A (I, J) coefficients --statements 
12 to 14- -with the appropriate terms from the desired expression. 

The program, as listed, accommodates 10 terms, but is easily altered 
to handle any desired number of terms. 

This program, as listed, requires an NBS fluid property data 
deck (data processing cards) for hydrogen or nitrogen fluids. 
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Computer Program Listing 


* JOB *2750453 *HORD > 2 
• FTN*L»X,R 

PROGRAM PO VO TO 

DIMENSION P ( 200 ) ,V(20 0) ,T<200) ,L(200) » PO ( 50 ) , VO ( 50 ) » TO ( 50 ) *K I V ( 50 ) 
1 » METPO ( 50 ) ,METT0(50) »METVO(50) 

COMMON A ( 200,10) ,3(20G),NP,NF,C(10*10),Y(10),D(10),R(75) 

DOUBLE PRECISION A»B»C»D»R»Y 
REAL KIV.METP0*METTO,METVO 

19 READ 72*1 CASE 

I F ( EOF >60)92*13 
13 K.TR= I = 1 
NF = 8 
SS = 0. 

10 READ 1 *L ( I) *T ( I ) * V ( I ) »P ( I ) 

I F ( EOF » 60 ) 22,11 

11 B ( I ) =P ( I ) 

SS=B ( I )*B( I ) +SS 

20 GO TO(12, 15,16, 17,18) ,ICASE 


12 A( I * 

1 ) =V ( I ) 

Ad* 

2)=T( I > 

A(I, 

3 ) =V ( I ) *V ( I) 

A(I * 

4 ) =T (I ) *T ( I) 

Ad* 

5)=T( I )*T( I)*T( I ) 

Ad* 

6 ) =V ( I )**4*T( I ) **2 

A ( I* 

7 ) =V ( I ) **4 

A ( I * 

8 ) = 1 • 00 

GO TO (14*44*55) *KTR 


15 A(I 


1 )=T(I ) 

A(I 


2 ) =T (I ) *V ( I ) 

Ad 


3 ) =T (I ) *T (I ) 

Ad 


4 ) =V ( I ) *V( I ) *T ( I) 

Ad 


5 ) = T ( I ) *T( I >*T ( I ) 

Ad 


6 ) =V ( I )**4*T( I ) **2 

Ad 


7 ) =V ( I ) **4 

Ad 


8 ) = 1 • 00 

GO 

T0( 14,44*55) »KTR 


16 Ad 

9 

1 )=T( I ) 

A ( I 

* 

2 ) =V (I > 

Ad 

9 

3 )=T(I )*V(I ) 

A ( I 

9 

4 ) =T ( I )*T< I ) 

Ad 

9 

5 ) =V ( I ) *V ( I ) 

Ad 

9 

6)=(T( I )*V( I ) )**2 

Ad 

9 

7 ) =V ( I ) **4 

Ad 

9 

8)=1.00 

GO 

T0( 14,44,55) ,KTR 


17 Ad 

9 

1 )=T(I ) 

Ad 

9 

2 ) =V( I ) 

Ad 

9 

3)=T( I )*V( I ) 

Ad 

9 

4 ) =V ( I )*T( I )*T( I ) 

Ad 

9 

5 ) =T ( I )*V( I )*V( I ) 

Ad, 

6)=(T( I )*V( I ) ) **2 

A ( I , 

7 ) =V ( I)**4 

A(I , 

8 1 = 1.00 

GO 

TO ( 14,44*55) »KTR 
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18 All* 1 )=V< I ) 

A ( I » 2)=T< I ) 

A ( I ♦ 3)=V< I >*V( I )*T( I > 

AC I , 4 ) =V ( I ) *V ( I ) 

A ( I ♦ 5)=T( I )*T< I ) 

A ( I > 6) = (T( I )*V( I ) ) **2 
A( I * 7 ) =V ( I ) **4 
A ( I » 8 ) =1 • 00 
GO TO ( 14 *44 > 55 ) » KTR 
14 1=1+1 

GO TO 10 
22 NP = I -1 
KTR = 2 

CALL BJORCK 
DO 24 I = 1 * 10 
S5 = SS-Y ( I )*Y( I ) 

FM=NP— I 

D ( I ) =SQRT ( SS/FN ) 

R(I)=D(I)*(C<( < 1. 870027491 /FN+1. 062 820 556 ) /FN+2. 63 502 131 >/FN 
1 +2.800533017) /FN+2. 377926828) /FN+1. 95967071 1 ) 

24 CONTINUE 

PRINT 2»(I*Y(I)*R(I)*D(I)»I=1»10) 

30 READ 6»NF*TZ*TM»DT *VZ»VM»DV*KK 
I F ( EOF , 60 ) 9,33 
33 PRINT 3 » ( C ( NF » J ) ,J=1,NF) 

PRINT 4 , N F 
DO 47 1=1 »NP 

GO TOC 12,15,16,17,18 ), ICASE 

44 PC=0. 

DO 45 J=1»NF 

45 PC=PC+A( I , J)*C(NF,J) 

D IFF = P ( I)-PC 

47 PRINT 5,L( I) *T( I) »V( I) ,P( I ) »PC»DIFF 
KTR = 3 
I =NP+1 
LL= 1 
TCI ) =TZ 

50 PRINT 7 , T ( I ) , NF 
V( I )=VZ 

GO TOC 12,15,16,17,18 ) »ICASE 
55 PC=0. 

DO 57 J=1»NF 
57 PC=PC+A( I ,J)*C(NF,J) 

PRINT 8, VC I > ,PC 
VOCLL )=V( I ) 

PO(LL)=PC 
TO ( LL ) =T ( I) 

LL=LL+1 

vc i )=vc d+dv 

IFC VC I ) .LE.VM) GO TO 20 

TCI) =T ( I)+DT 

IFCTC D.LE.TM) 50,30 

1 FORMATC A8,F12.0,2F10.0) 

2 FORMAT ( * 1 THE NUMBER OF TERMS, INTERNAL COEFFICIENTS, THEIR ESTIMAT 
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1 ED ERROR* AND THE STANDARD DEVIATION FOR THAT MANY TERMS ARE AS FO 
2LLOWS9*/*0*/{*OA(*I2»*) =* E17.10** +/-*E9.2,*, STD. DEV. =*E9.2)) 

3 FORMAT ( *1CALCULATED PO is THE SUM OF THE FOLLOWING TERMS9*/*0*/ 

1 /*0 PO =*E1 7. 10 » 3H*A ( I > 1) E17.10,8H*A( I ♦ 2) E17. 10»8H*A( I * 3) 

2 /*0 *E17. 10 * 8H*A ( I * 4) El 7 . 10 . 8H*A ( I * 5) E 1 7. 10 , 8H*A ( I , 6) 

3 /*0 *E17.1C,8H*A( I » 7) E 1 7 . 10 * 8H*A ( I , 8) E17. 10»8H*A( I » 9) 

4 /*0 *E17.10,3H*A( I ,10) E 1 7 . 10 , 8H*A ( I , 1 1 ) E17. 10 ,8H*A( I , 12 ) 

5 /*0 *E17.10,8H*A( I *13) E 1 7. 10 , 8H*A ( I , 14 ) E17. 10 , 8H*A ( I , 15 ) 

6 /*0 *E17.10» 8H*A (1*16) E17. 10 ,8H*A( I *17 ) E 17. 10 , 8H*A ( I , 1 8 ) 

7 /*0 *E17. 10 »8H*A ( 1*19) El 7 . 10 ♦ 3H*A ( I * 20 ) E 17. 10 , 8H*A ( I , 2 1 ) 

8 /*0 *E17. 10 ,8H*A ( I *22 ) E 1 7 . 1 0 * 8H* A ( I » 23 ) E17.10»8H*A( I ,24) 

9 /*C *E17.10»SH#A( I *25 > ) 

4 FORMAT ( *1*12 ** FUNCTIONS GIVES THE FOLLOWING FIT9*/*0*/ 

1 * LAEEL T V P P CALC DIFF*/* *) 

5 FORMAT ( * *A4,F8.2»F7.1*F8.1»F8.1>F6.1) 

6 FORMAT! 110, 6F10.0, 110) 

7 FORMAT ( *1THE*F6. 2 ** DEGREE R. ISOTHERM USING*I3»* FUNCTIONS IS9* 

1 /*0*/* V PO *) 

8 FORMAT (*0*2F8.1) 

9 LL=LL-1 

60 IF(KK.EQ.1)90,61 

C THIS BRANCHES TO KIV CALC IF KK NOT EQUAL TO 1 

61 DO 62 J= 1 »LL 

METPO( J)=0.6894757*PO( J) 

METTO( J)=(5./9. ) *T0 ( J ) 

MET VO ( J ) =0. 3048*VO ( J ) 

ATMPO=PO( J) / 14. 6959494 
PSAT =VPRESS ( MET TO ( J ) ) 

RH0LIQ=FIND DH( ATMPO,METTO( J) ) 

I F ( RHOL IQ.GT.O.O) GO TO 63 
K I V ( J ) =0 » 0 
GO TO 62 

63 KIV(J)=(PO( J)-14.6959494*PSAT)*144.*2.*32.1725/(RHOLlQ*VO( J)**2) 

62 continue 
WRI TE ( 61 ♦ 70 ) 

WRITE ( 62 *70 ) 

WRITE (61*71 ) (TO(J) ,VO(J) ,PO(J) ,KIV(J) ,METPO ( J ) ,METVO ( J ) ,MeTTO( J) , 
1 J=1 »LL ) 

WRITE (62*71 ) (TO( J) *VO( J) ,PO( J) ,KIV( J) ,METPO ( J ) ,METVO ( J ) ,MeTTO< J) * 
1 J=1 ,LL ) 

70 FORMAT ( * 1* , 4X,*TEMPERATURE~C0MPENSATED DESINENT DATA ( HYDROGEN90 . 4 
120-INt 1.067-CM) OGIVE)* /// 7X , *TO* , 8X »*VO* » 8X ,*PO* , 10X » *K I V* , 10 
2X,*PO*,8X,*VO*,8X,*TO*/ 6X,*DEG R*,4X»*FT/SEC*»5X»*PSIA*,20X,*N/CM 
3/CM*,4X,*M/SEC*,5X,*DEG <*//) 

71 FORMAT ( 6 ( IX , F 1 0 . 2 , F9 . 1 , F 1 1 . 2 » FI 2 . 2 , FI 3 . 2 ,F9. 1 ,F 1 1 . 2 , / ) / ) 

72 FORMAT (II) 

90 GO TO 19 
92 CALL EXIT 

END 

FUNCTION FIND DH(PATM»TDK) 

CALCULATES THERMOFUNCTIONS OF H2 ACCORDING TO MONO 94 
INPUT , P IN ATM, T IN DEG K , OUTPUT UNITS AS IN MONO 94 
EXCEPT THAT DENSITY IN MOLES/CC IS RETURNED, NOT VOLUME, AND 
INTERNAL ENERGY HAS TO BE CALC. FROM H AND P-V. 
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INPUT LIMITATIONS, P LESS THAN 350,1, T BETWEEN 13.803 AND 100.0 
SUBROUTINES REQUIRED, THERMO, CASEI ,CASEI I ,CASEA, IS0THM»C2l4TApE 
DIRECT ANY QUESTIONS TO HANS RODER 

COMMON /C214TAPE/CT(43? 15) ,CR<90,5 ) >RR(90) ,TT (43) ,CW( 100) »SW( 100) * 

1 HW( 100 ) , PST ( 2 8,8) ,TTT( 66 ) ,-CVSPt 66) ,C ( 2 ,6 ) ,PP ( 90 ) ,DPDT( 90) , 

2 DPDT2 (90) ,B1,B2,B3,B4,CA,CB»CC,CD,CE,DPDR,DPDR2,DPDR3,PI,MI,I I , 

3 JJ» I STOP, istopi»jstop,cpdtdr»cpdtdi»p*t»rl,rv,h,5»cv»cp»svel» 

4 CSL0PE»P1ST (28,7) ,P2ST (28,7) ,P3ST (28,6) 

P=PATM 

T = TDK 
M 1 = 1 

CALL CASE A 

FIND DH=RL*125. 83728 

RETURN 

END 

SUBROUTINE BJORCK 

COMMON A ( 200,10) *B ( 200 ) ,NP ,NF ,C ( 10 , 10 ) , Y ( 10 ) , D( 10 ) ,R ( 75 ) 

DOUBLE PRECISION A,3,C,D*R»Y 
D ( 1 ) =Y 1 1 ) =C »D 
DO 1 1=1, NP 

D ( 1 ) = A ( I » 1 ) *A ( I , 1 ) +D ( 1 ) 

1 Y( 1 )=A( 1 ,1>*B( I )+Y( 1 ) 

Y( 1 )=Y( 1 ) /D< 1 ) 

I R = 0 

DO 5 K = 2 »NF 
DO 3 J=K»NF 
I R= I R+l 
R( IR)=O.D 
DO 2 1=1, NP 

2 R(IR)=A(I,K-1)*A(I,J)+R(IR) 

R(IR)=R(IR)/D(K-1) 

DO 3 1=1, NP 

3 A( I ,J)=A( I, J)-A( I ,K-1)*R( IR) 

D ( K. ) =Y ( K.) =0»D 

DO 4 1=1, NP 

B( I)=B( I)-A( I »K~1 ) *Y ( K— 1 ) 

Y(K)=A( I » K ) *B ( I ) +Y ( K ) 

4 D ( !< ) = A ( I *K)*A ( I ,K)+D(K) 

5 Y ( K ) =Y ( K) /D l K ) 

IRS=-NF 

DO 7 K= 1 , NF 
IR=IRS=I RS+NF-K+1 
DO 7 JJ=1,K 
J=K-JJ+1 
C ( K , J ) =Y( J ) 

DO 6 I =2 , JJ 

C(K,J)=C(K,J)-C(K,K-I+2)*R( IR) 

6 I R= I R— 1 

7 IR=IR-NF+K 
DO 8 1=1, NF 

8 Y( I )=Y( I )*D( I ) **0 • 5 
RETURN 

END 

SCOPE 


-7111 
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